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ABSTRACT 

We examine the changes in the properties of galactic bulges and discs with environ- 
ment for a volume-limited sample of 12500 nearby galaxies from SDSS. We focus on 
galaxies with classical bulges. Classical bulges seem to have the same formation his- 
tory as ellipticals of the same mass, and we test if environment determines whether 
or not a classical bulge possesses a disc. Using the projected fifth nearest neighbour 
density as a measure of local environment, we look for correlations with environment 
at fixed bulge stellar mass. In groups with fewer than 20 members, we find no evidence 
for changes in disc morphology with local density. At fixed bulge mass, disc mass and 
disc scale length are independent of local density. However, disc colour does increase 
(A(g — r) ^ 0.05 mag) as a function of local density in relatively poor groups. There- 
fore, the colour-density relation for classical bulge-|-disc galaxies in the field and in 
poor groups is due solely to changes in disc colour with density. In contrast, we find no 
correlations between disc colour and local density for classical bulge-|-disc galaxies in 
large, relaxed groups and clusters. However, there is a weak correlation between disc 
mass and group crossing time, suggesting morphological transformation takes places 
in rich groups. Our results add to the evidence that star formation is quenched in 
group environments, instead of clusters, and that star formation quenching and mor- 
phological transformation are separate processes. Overall, we show that environment 
has two effects on galactic discs: relatively low density environments can quench star 
formation in discs, while processes occurring in higher density environments contribute 
to the morphological transformation from disc-dominated systems to bulge-dominated 
systems. 
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1 INTRODUCTION 

Galaxy morphology, stellar mass, star formation rate, 
and projected number density are all known to cor- 
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relati vely low density regions (e.g. Mulchaev fc Zabludofll 
[l99i). For galaxies which follow the general pattern, it is 
unclear which of the correlations mentioned above are the 
result of physical processes, and which, if any, are simply 
consequences of other correlations. Morphology, star for- 



lar mass (e.g. HamiltonI 19881: Brinchmann & EUisI 2000l: 
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are disc-dominated, star-forming (blue), and reside in low 
density regions. There are, of course, exceptions: passive 
discs make up a significant fraction of the r ed sequence 



(e.g. iBernardi etHI l2003al : iMaller et all |2009|) blue eUip- 
ticals are still forming stars ( Schawinski et al] |2009|) , and 



these properties are partially due to their correlations 
with stellar mass. Yet, at fixed stellar mass, stud- 
ies have found correlations between density and mor- 
phology (^B amford ct al. 200^), density and stellar age 
(Tho mas et al.i. 2005: Co oper et al.ll2O10l), density and colou r 
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((Balogh et al.ir2004. : ,Skibba et al.ll2"009l: ICibinel et al.ll2012l) 
and density and star formatio n rate (jKauffmann et al .I I2OO4I : 
IChristlein fc Zabludofj I2OO5I ) . Furthermore, at fixed lumi- 
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nosity, neither blue nor red galaxy colours seem to depend on 
density, i.e. blue galaxies do not get redder as a function of 
density, only their num ber fraction decreases (|Balogh et al.l 
I2OO4I : iHogg et aLlbOOal l . Furthermore, correlations with mor- 
phology aB^_d«nsity seem to disappear for high m ass galax- 
ies (e.g. lTasca et ai.ll2009l : iGriitzbauch et al.li201ll ). 

The local correlations between galaxy properties and 
density extend to higher redshift; the morphology-density 

and colour-density relations are in pl ace by 2: » 1 

llDressler et al.lll997l : iPostman et"al]|2005l : iTreu et al.ir2003l : 
ISmith et al.ll2005l ). but the relations do evolve with redshift. 
The fraction of blu e galaxies in high density regions increases 
with redshift (e.g. iButcher fc Oemlei|[l978l '). while the frac- 
tion of SOs and red discs galaxies decreases with increa sing 
redshift CSniith ot al. 2005; Moran ct al. 2007; Bundy et~al] 
|2010I :Ib ruce et al.ll2012l ). In addition, the hierarchical growth 
of structure implies that galaxies generally move from low 
density regions to high density regions as a function of time. 
Therefore, any environmental effects and trends will be more 
pronounced today than in the past (jTasca et al.| |2009i). 

From these observations, a general outline of galaxy 
evolution has been developed. At early times, galaxies are 
blue, intensely star-forming, disc systems. As galaxies be- 
come more massive, star formation is quenched (by internal 
feedback mechanisms), creating a population of re d, mas- 
sive galaxies. This mass quenching (|Peng et al.ll2010l ) is com- 
pounded by environmental effects. At fixed mass, galaxies in 
higher density regions become red and bulge-dominated ear- 
lier, creating the colour-density and morphology-density re- 
lations. Star formation quenching is thought to occur before 
morphological transformation, which leads to an increase 
in the fraction of SOs at intermediate densities (e.g. 'DressW 
1980; Mcintosh ot al. 2004; Cooper ct al. 2006; Moran ct al. 
20071 : iBundv et al.ll2006l ). In order to separate the effects of 
environment from the effects of stellar mass, we examine 
correlations between galaxy morphologies, colours, and lo- 
cal density at fixed stellar mass. 

The physical processes responsible for the environment- 
driven transformations are unknown, although there are 
many candidates. Processes can be divided into those which 
truncate star formation, and t hose which also cause m or- 
phological transformations (see lBoselli fc Gavazzi|[2006l . for 
a review of these processes). R am-pressure strippin g of ISM 
from galaxies entering clusters (|Gunn fc Gottlll972[) and the 
removal of hot halo ga s (strangulation) ( Larson et al.lll980l : 
iBalogh fc Morris! I2OOOI ) both act to truncate star formation 
in discs and can transform spiral discs into SOs, but do not 
drastically alter a galaxy's stella r disc. Tidal stripping by the 
cluster potential (jMerrittlll98i ) and hig h speed encounters 
with other cluster galaxies (harassment) l|Moore et al.lll99^ . 
Il998l . Il999l ) both act to transform disc-dominated galaxies 
into bulge-dominated galaxies. All these processes act on 
a variety of timescales and require different minimum lo- 
cal densities in order to be effective. It is likely that more 
than one process is responsible for the morphology-density 
relation. 

In this work, we examine low redshift galaxies with both 
a bulge and disc, and study the correlations of the separate 
bulge and disc properties with local density. These galaxies, 
which include SOs, may represent a transition from disc- 
dominated to bulge-dominated, and the environmental pro- 
cesses enumerated above should have observable effects on 



the discs and possibly the bulges of these transitional galax- 
ies. For the bulge and disc properties, we u se the bulge+disc 
decom positions from our earlier work (|Lackner fc GunrJ 
|2012| . hereafter L12). L12 presents bulge-|-disc decomposi- 
tions for nearly 72,000 low redshift (0.002 < z < 0.05) 
galaxies from the Sloan Digital Sky Survey (SDSS). The 
galaxies we use for in this work are a luminous subsample 
from L12. 

We focus on galaxies which host classical bulges. 
These bulges have properties and, presumably, formation 
histories identical to elliptical galaxies of the same mass. 
Classical bulges are concent r ated, p ressure-supported 



syste ms (|Falc6n-Barroso et al.l |2002|: Mac Arthur et all 
2008h. with old st ellar populations ("Peletier et al. 19991 : 



Moorthv fc Hohzman 2006; MacArthur ct al. 20T^, but see 



Gadotti 



LV &: ±lc 
il l2009l ') 



In L12, we model the light profiles of both 
classical bulges and elliptic als using a de Vaucouleurs profile 
(but see ICaon et al.|[l993l '). We show in L12 that classical 
bulges and ellipti cals follow the same size-density relation 
(Kormcnd^ ll977l ). Taken together, these studies support 
the assertion that classical bulges and ellipticals of the same 
stellar mass are indistinguishable. Therefore, by exploring 
how local density affects the discs around classical bulges, 
we can determine if changes in environment correspond 
to a transition from classical bulges with discs to disc-less 
elliptical galaxies. 

To date, there have been a handful of studies which 
explore the effe cts of environment on bul ges and discs 
separa tely. Both iMcIntosh et al.l (|2004l 'l and [Hudson et af] 
ll2010t) show that disc colour is a function of cluster radius. 
iMcIntosh et~all (|2004 ) also show that the amount of sub- 
structure in discs declines with increasing density, further 
demonstrating that star formation is quenched in dense en- 
vironments. The sample we use in this work is considerably 
larger than the samples in previous studies, and it covers 
the entire spectrum of local densities, not just rich clusters 
and the field. Furthermore, our large sample can easily be 
divided into subsamples of constant bulge mass, eliminat- 
ing trends with stellar mass and environment, and still yield 
statistically significant results. 

Below, we use our sample of classical bulge-f disc galax- 
ies to determine whether the colour-density relation for 
these galaxies is due to changes in bulge or disc colour, 
indicating star formation truncation, or changes in bulge- 
to-total ratio, indicating morphological transformation. We 
find that the colour-density relation for these galaxies at 
fixed bulge mass is due entirely to changes in disc colour, not 
changes in disc mass or size. Next, we divide the sample into 
rich and poor groups to determine if the trends in disc prop- 
erties with environment depend on group size or halo mass. 
We find that while disc colour is a function of local density 
in relatively poor groups, disc colour is independent of local 
density in larger, relaxed groups and clusters. However, Disc 
mass decreases slightly with increasing local density in large 
groups and clusters, while disc mass is independent of local 
density for galaxies in the field and in poor groups. From 
these results, we conclude that environment-driven star for- 
mation quenching occurs in relatively low density environ- 
ments, while structural changes to discs only occur in higher 
density environments. 

In order to perform the studies detailed above, we re- 
quire a robust sample of classical bulge-|-disc galaxies. It 
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is especially important to distinguish classical bulge+disc 
galaxies from disc-less ellipticals. We present a probabilistic 
method for separating classical bulge hosts and ellipticals us- 
ing bulge-fdisc decompositions in Section [2. 1.1 1 The remain- 
der of Section [2] details the environment metrics employed. 
In order to determine group membership and local density, 
we use an updated group catalogue from A. Berlind, (priv. 
co mm.), which is base d on the group catalogues presented 
in lBerlind et al.l (|2006l ). Section [3] presents the morphology- 
density relation and colour-density relation for our sample. 
The results of this section help to confirm the assignment of 
galaxy morphologies in Section 12.1.11 Sections |4] and [5] are 
devoted to the correlations between environment and bulge 
and disc properties. Section |4] focuses on correlations with 
projected local density, while Section [5] discusses classical 
bulge-|-disc galaxies in relatively rich groups. 

Throughout this paper we use the ACDM cosmology: 
= 0.3, Ho = 70 km s~^Mpc~\ and D,x = 0.7. 



2 SAMPLE 

The sample used in the work consists of a luminous subsam- 
ple of galaxies from our earlier bulge-|-disc decompositions 
matched to an updated group catalogue from A. Berlind 
(priv. com m.). This group cata logue employs the method 
presented in iBerlind et alj(|2006t ). but uses data from SDSS 
data release 7. The sample contains 29781 galaxies and is 
complete to an absolute "'^r— band magnitude of —19.77 and 
covers the redshift range 0.02 ^ z < 0.05. Below, we describe 
the bulge-|-disc decompositions and the environmental infor- 
mation extracted from the group catalogue. 

2.1 Bulge-Disc Decompositions 

For the bulge and disc properties, we use the results of our 
earlier work (L12). L12 presents bulge-fdisc decompositions 
for 72, 000 galaxies from SDSS data release 8 (the data is 
the same as that in data release 7, but the reductions have 
been improved). The galaxies have redshifts between 0.003 
and 0.05. All of the galaxies are in the SDSS spectroscopic 
sample, which implies a limiting magnitude of < 17.77. 
Two dimensional bulge-disc decompositions are performed 
for the r— band images. The results are then linearly scaled 
to fit the galaxy images in the u, g, i, and z bands, yielding 
colours for the bulge and disc components. Because we only 
linearly scale the fits in each band, our bulge-fdisc mod- 
els do not take into account colour gradients within each 
component. Each galaxy is fit with 5 different models: a de 
Vaucouleurs bulge and an exponential disc (rib ~ 4 B+D), 
an exponential bulge and exponential disc (rit = 1 B+D), 
a single de Vaucouleurs profile, a single exponential profile, 
and a single Sersic profile. 

The two bulge-fdisc models allow us to fit both 
elliptical-like, pressure-supported classical bulges (de Vau- 
couleurs profile) and disc-like, rotationally-supported 
pseudo-bulges (expon ential profile) (|Kormendvlll977l . 1 19931 : 
iFisher fc Drorvl [20081 ) . Pseudo-bulges are thought to arise 
from secular p rocesses within discs, such as bar-driven in- 
stabi li ties (e.g. [Ko rmendv fc Ke nnicutt|[2004l : lAthanassoulal 
I2OO5I : IWeinzirl et al. ^09), and, as such, have very dif- 
ferent formation histories than classical bulges (but see 



lElmegreen et a l."2009'). Often, pseudo- bulges arc still form- 
ing stars today (Kormendv fc Kennicutt,200 4: Fisher 200Q) . 
Since pseudo-bulges are a disc phenomenon, we do not in- 
clude them in our sample of bulge-fdisc galaxies. 

For each galaxy, we tabulate the bulge and disc mag- 
nitudes and colours in all 5 S DSS ban ds. These values are 
Galactic extinction corrected (|Schlegel ct al . 1998) and k- 
corrected to z = usin g the IDL package kcorrect v4_2 
(IBlanton fc Roweisll2007l ). In addition, we correct the colours 
and magnitudes of gala xies with disc s for in trinsic extinction 
using corrections from iMaller et all (|2009l ) and L12. These 
corrections remove trends in colours with disc inclination, 
but they do not correct for extinction due to dust in face-on 
discs. Finally, we calculate the stellar mas s es for the bulge 
and disc using the relation from lBell et al.l (|2003l ): 

logM/M© =- 0.22 + 0.66(5-4) -0.15 (1) 
- 0.4(M^ - M©,,. + 1.32) , 

where 0.15 accounts for the difference betwee n the diet 
Salpetcr initial mass function (IMF ) use d by iBell et al] 
( 2003) and the Kroupa IMF ( Kroupa 20021 ) we employ. The 
colours and magnitudes used for the stellar mass are not cor- 
rected for intrinsic extinction , in kee ping with the derivation 
of the relation in iBell et all (|2003h . Because the mass-to- 
light ratio is a convex function of galaxy colour, the sum 
of the masses of the bulge and disc is always slightly larger 
than the mass measured using the total galaxy colour and 
magnitude. For most of the galaxies, this difference is small; 
the median Mtotai/(Mdisc + Mbuigc) is 0.9996, and for 95 per 
cent of the galaxies, this ratio is between 0.85 and 1.0. 



2.1.1 Classifying Galaxies 

Since we fit each galaxy with five different models, we re- 
quire a method for selecting the best-fitting model. In L12, 
we show that the values of the various model fits ar e indis- 
tingu ishable at the resolution of SDSS (but see .Simard et al.l 
I2OIII ). Instead, we develop an algorithm that relies on the 
sizes, shapes, and colours of the bulges and discs in order 
to select the best-fitting, physically-sensible model for each 
galaxy. In this work, we present a simplified classification 
algorithm that attempts to classify most of the galaxies in 
our sample and emphasises the distinction between classical 
bulge galaxies and ellipticals. Additionally, instead of assign- 
ing each galaxy a best-fitting model, we assign each galaxy a 
probability of being fit by each model (often, the probability 
is unity for one of the models). Although this does not al- 
low us to accurately classify a given galaxy, it does allow us 
to study the properties of a large sample of galaxies. When 
we examine properties of bulges and discs, we weight each 
galaxy by its probability of having a bulge and a disc. 

We separate galaxies into five different categories: 
bulge-less disc galaxies, disc-less ellipticals, classical 
bulge+disc galaxies, pseudo-bulge+disc galaxies, and un- 
classifiable galaxies. Our goal is to assemble a sample of 
galaxies which are accurately modelled by a classical bulge 
plus a disc. A brief outline of the classification is given be- 
low and summarised in Table \T\ Details can be found in 
Appendix [XI 

First, we identify bulge- less and disc-less galaxies. 
Bulge- less disc galaxies are defined to have B /T < 10 per 
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Table 1. Classification scheme for galaxies. The number of galax- 
ies refers to the number in the bright, grouped sample. For most 
of this work, we only use the classical bulge+disc galaxies. 



Type 



Num 



Classification 



elliptical 



3681 



either B/T > 0.9 or bulge and disc 
u — r > 2.22 and statistically sepa- 
rated from red B-l-D galaxies based 
on disc axis ratio (Fig. IA2|| 
classical B+D 12523 if bulge and disc u — r > 2.22, sta- 
tistically separated from ellipticals 
based on disc axis ratio (Fig. lA2ll : 
otherwise, selected to be quiescent 
based on Z)„(4000) to distinguish 
from pseudo-bulge hosts(Fig. IA3|| 
pseudo B+D 5361 bulge or disc ti — r ^ 2.22 and 
statistically separated from classi- 
cal bulge hosts based on D„(4000); 
pseudo-bulges are selected to be 
star-forming 
bulgeless disc 8505 B/T < 0.1 

unclassifiable 2857 none of the above categories, typi- 
cally blue, disc-like irregulars 



cent. Disc-less galaxies (ellipticals) are more difficult to iden- 
tify. The few galaxies with B /T > 90 per cent are considered 
ellipticals. These make up only 6 per cent of th e our sam- 
ple. As shown in L12 (see also I Allen et al.ll2006l ). elliptical 
galaxies are often best fit by a de Vaucouleurs component 
along with a low surface brightness exponential component. 
This "disc" is not a physical disc and has several possible 
origins, i.e. the outer halo of ellipticals, a Sersic index larger 
than four, and/or inadequate sky subtraction around bright 
galaxies in SDSS. These model "discs" make ellipticals in- 
distinguishable from face-on bulge-|-disc galaxies based on 
the 2-dimensional bulge-|-disc decomposition alone. 

However, the distributions of inclination angles for real 
and spurious discs will be different; the former will be ran- 
domly oriented, while the latter will be preferentially face- 
on. We use this fact to statistically separate ellipticals from 
face-on classical bulge host galaxies. Galaxies with a small 
measured disc axis ratio have a high probability of being 
a bulge-|-disc galaxy, while galaxies with a large disc axis 
ratio (face-on) might be either a bulge-|-disc galaxy or an 
elliptical with an exponential halo. For this statistical sep- 
aration, we examine galaxies with 0.1 < B/T < 0.9 and 
{u — r) > 2.22 for both the bulge and disc component. This 
means all the ellipticals we find will be red. Restricting our- 
selves to red galaxies will enhance the fraction of ellipticals 
relative to bulge-|-disc galaxies, making the two inclination 
angl e distributions easier to fit. Blue ellipticals are relatively 
rare l|Schawinski et al.ll2009l ). and, therefore, are a small con- 
tamination in our sample of classical bulge host galaxies. 

After setting aside ellipticals, we distinguish between 
galaxies with quiescent bulges and those with young, 
star-forming bulges based on the 4000A break strength 
{D„{4000)) measured by SDS^|. Based on results from mor- 
phology studies at higher resolution (e.g. iFisherl [20061 ) ■ we 
associate star-forming bulges with pseudo-bulges and quies- 



^ http:/ /www. mpa-garching.mpg.de/SDSS/DR7/raw_data. html 



cent bulges with classical bulges. As above, we assign each 
galaxy a probability of having a classical bulge or pseudo- 
bulge based on its D„(4000) (see Fig. IA3|l . Since pseudo- 
bulges are a disc phenomenon, we consider galaxies with 
pseudo-bulges to be bulge-less and exclude them from our 
analysis. Because 90 per cent of the star-forming bulges are 
less massive than the median classical, quiescent bulge in our 
sample, excluding pseudo-bulge hosts does not significantly 
affect correlations in galaxy properties with local density at 
fixed bulge mass. However, the separation of quiescent and 
star-forming bulges will eliminate any star-forming classical 
bulges from our sample. 

Finally, we exclude galaxies which are not well fit by any 
of the above models. These galaxies are modelled by a single 
Sersic profile. Unclassifiable galaxies make up less than 10% 
of our sample. Seventy-five percent of unclassifiable galaxies 
have a Sersic index less than 2.3 and the same fraction lie in 
the blue cloud (u — r < 2.22). These galaxies are probably 
disc-like irregulars, which are unlikely to have a well-defined 
bulge and disc. The remaining 25 per cent of unclassifiable 
galaxies are mostly merger remnants, starbursts, and other 
complicated morphologies. None the less, because the ma- 
jority of unclassifiable galaxies exhibit disc-like properties, 
we group them with other bulge-less galaxies. 

Below, we use a bright (Mo.i^ < —19.77) subsample of 
29781 galaxies from the L12 sample. This sample contains 
12523 classical bulge-|-disc galaxies, 3681 ellipticals, 2214 
pseudo-bulge galaxies, 8505 bulge-less galaxies, and 2857 
unclassifiable galaxies. These numbers are inexact, since we 
only assign each galaxy a probability of being a certain type. 



2.2 Group Catalogues 

We study the environmental properties of bulges and discs 
using the group catalogue from A. Berlind ( priv. comm.). 
This catalogue is built using the methods in iBerlind et al.l 
(2006), but is based on SDSS data release 7 instead of data 
release 4. The group catalogue allows us to relate galaxy 
properties to their host group (and dark matter halo) prop- 
erties as well as to study the properties of galaxies as a 
function of intra-group environment. We select this group 
catalogue from the many group catalogues available for 
SDSS data because it extends to low redshift, and overlaps 
significantly with our sample of bulge-|-disc galaxies. The 
group catalogue is volume-limited and includes all galax- 
ies with absolute magnitudes Mo.i^ ^ —19.77 in the red- 
shift range 0.02 < z < 0.067. We include isolated galaxies 
in the group catalogue (groups of richness one). The cata- 
logue is created using a friends -of-friends algorithm to de- 
termine group membership (see lBerlind et aL I l2006l . for de- 
tails). Two galaxies are linked if the projected and trans- 



verse distances between them are smaller than bin 



— — 1 /3 

feynj , respectively, where 6|| x are the linking lengths and 
fig is the average galaxy number density in the sample. Here, 
fig = 5.275 X 10"^ Mpc-^ 6|| = 0.14, and 6x = 0.75, which 
corresponds to physical linking lengths of 0.8 Mpc in the 
transverse direction and 300 km s~^ along the line of sight. 

In the group catalogue, there are 90893 galaxies, of 
which 29781 have a bulge-|-disc model fit from L12. This sub- 
sample includes 83 per cent of the galaxies below z — 0.05 
with spectroscopic redshifts in the group catalogues. Galax- 



-1/3 



and 
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S5 [Mpc-2] 

Figure 1. Top panel: The distr ibution of galaxy gro up sizes for 
the total group catalogue from iBerlind et alj |200g) (blue solid 
line) and the subsample with bulge-disc decompositions (black 
dashed line). Lower panel: The distributions of S5 for the same 
samples. These two distributions are indistinguishable at the 2a 
level. 



ies without spectroscopic redshifts (due to fibre collisions in 
the SDSS spectrograph) make up 4 per cent of the sample 
below z — 0.05, and are therefore a small omission. The re- 
maining 6124 galaxies are missing from the L12 bulge+disc 
sample because cuts in the galaxy axis ratio (1289), model 
fits with surface brightness consistent with zero (2830), and 
galaxies which did not make our quality cuts due to prob- 
lems with deblending and cosmic ray removal. 

Despite the missing galaxies, the bulge+disc matched 
sample is a representative subsample of the group catalogue 
for z ^ 0.05. This is demonstrated by Fig. [1] which shows 
the distributions of group sizes and projected fifth nearest 
neighbour densities (E5, see i|2.3l for explanation) for the full 
group catalogue and for the subsample which overlaps with 
the bulge-|-disc decomposition sample. The bulge+disc sam- 
ple follows essentially the same distributions in E5 and A'^gai 
as the full group catalogue. A Kolmogorov-Smirnov (KS) 
test shows that the distributions of E5 (lower panel) are 
indistinguishable. The distributions of group sizes are not 
identical; a KS test yields a probability of 3 x 10~^ that the 
distributions are the same. The bulge+disc matched sample 
is missing galaxies in mid-sized groups, but the number of 
missing galaxies is small. Furthermore, for the larger groups, 
the missing galaxies are not a function of position in the 
group or of local density. This is not surprising, since the 
number of galaxies missing a spectrum due to fibre colli- 
sions is a small fraction (4 to 6 per cent) of the sample. If 
fibre collisions were more prevalent, we would expect to be 
missing more galaxies in high density regions than in low 
density regions. 

We also check that matching the L12 sample to the 
group catalogue does not change the distributions of galaxy 
morphology as a function of magnitude. The top panel of 
Fig.[2]shows the absolute magnitude distribution of the total 




-20 -21 -22 -23 



Figure 2. Top panel: The distribution of galaxy absolute magni- 
tude for the galaxies from L12 with Mr < —19.77. The different 
colours represent the different categories of galaxies. The number 
of galaxies in each category is the sum of the probabilities calcu- 
lated in i|2.1.1l Bottom panel: The same distribution for galaxies 
from L12 matched to the group catalogue, which has a limiting 
magnitude of Mo.i^ = —19.77. 



L12 sample, down to a magnitude of Mr — —19.77. The 
lower panel shows the same for the galaxies from the group 
catalogue. Although there are small differences at the faint 
end, the distributions of galaxy morphology as a function 
of magnitude are essentially the same for the both samples. 
The fraction of galaxies of each type in the two samples 
is the same to within 450 galaxies, or 1.5 per cent of the 
grouped sample. 

In this work, we address the effects of environment on 
galaxy properties at fixed bulge mass. Therefore, although 
it is important that certain types of galaxies are not system- 
atically excluded as a function of environment, we are not 
concerned with the overall completeness of the sample. We 
are only concerned that the galaxies in the different environ- 
ments are a representative sample, which is demonstrated by 
Figs. [I and H 



2.3 Measuring environment 

In order to examine the correlations between galaxy prop- 
erties and environment, we use several measures of envi- 
ronment, including both group halo properties and more 
local measures of environment. Directly from the group 
catalogue, we obtain the group richness, Afgai, the group 
line-of-sight velocity dispersion, and the total stellar mass 
in galaxies brighter than Mo.i,, — —19.77. Using the rela- 
tion between tot a l grou p stellar mass and halo mass from 
iLeauthaud et al.l (|2012l ). we calculate the group dark mat- 
ter mass, M200, defined as the mass enclosed in a region 
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200 times denser than the critical densitj0. We do not take 
into account the difTerence in stell ar mass completeness be - 
tween our sample and that used in lLeauthaud et al] l|2012l ). 
nor do we correct the stellar masses for contamination from 
non- group galaxies. However, we expect the corrections to 
the stellar masses to be small fsee lLeauthaud et alj|201^ 1. 

We can compute the line-of-sight velocity dispersion (cr) 
for a halo o f a giv en mass using the Mhaio — o" relation from 
lYang et alj (|2007l ). For groups with more than 10 galaxies, 
the a obtained from the dark matter halo mass is a factor 
of 1.4 larger than the a measured directly from the galax- 
ies. Half of this discrepancy is due to the small value for the 
line-of-sight linking length used to build the group catalogue. 
The small fey biases the measured velocity dispersion down 
by ~ 20 per cent (Borlind ot al. 2006). In the following anal- 
ysis, we do not make extensive use of the dark matter halo 
mass. We do use the group velocity dispersion, but since we 
are only interested in making comparisons between different 
environments, the absolute values of a are not relevant. 

The group catalogue also contains the the projected dis- 
tance, Rp, from each galaxy to its host group centre (the 
number- weighted mean angular position). We use this dis- 
tance, along with the velocity dispersion of the group to 
define the crossing time tcross ~ Rp/o', where a is the the 
group velocity dispersion, measure d directly from the galaxy 
redshifts (see iBerlind et al]|2006l ). Clearly, tcross is only a 
sensible measure of environment for relatively large groups 
and clusters, which have a definite centre. While the number- 
weighted (or mass- weighted) centre of a group always is well- 
defined, the FoF grouping algorithm does not guarantee it is 
physically meaningful. This is especially true for non-relaxed 
groups. We will address this problem in our investigation of 
intra-group trends by limiting our sample to galaxies from 
relaxed groups ( H5.ip . 

Finally, in addition to the halo properties, we measure 
the local surface density around each galaxy. We use the pro- 
jected fifth nearest neighbour density (E5 [Mpc~^]). Neigh- 
bours are selected from the volume-limited group catalogue 
in a redshift slice which has a width equal to the velocity 
dispersion of a galaxy's host group. The minimum width is 
300 km s~^, the line-of-sight linking length used in the group 
catalogue. For galaxies in large clusters, E5 is a measure of 
their immediate vici nity, not the underlyin g large-scale dark 
matter density field (|Muldrew et al.ll2012l ). 

In the sections below, we rely mostly on E5 as a metric 
for environment. We will demonstrate that our results are 
essentially unchanged if group richness, A'gai, is substituted 
for E5. In i]5.1l we explore trends in galaxy properties within 
massive, relaxed groups and use the crossing time, tcross, as 
a measure of galaxy environment. 



3 WHOLE GALAXY PROPERTIES WITH E5 

Before examining trends in bulge and disc properties 
with density, we confirm trends in whole galaxy prop- 
erties with density. The low redshift morphology-density 
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Figure 3. The distribution of morphological types as a function 
of local density, S5. Classical B-l-D galaxies are the nj, = 4 B-l-D 
galaxies; pseudo- (star-forming) bulge galaxies are the = 1 
B-(-D galaxies. The unclassifiable galaxies are Sersic galaxies. The 
errorbars are the Poisson errors in the number of galaxies in each 
bin. 



and colour-density relations have been m easured numer 
ous times using various data sets (e.g. 



Goto et al.ll2003l:lBalo^Ii et al.ll2004 



Dressier 



De Propris et al 



1981 



200' 



iLeauthaud et al. report halo masses as 500 times the 

critical densit y. We convert to M2 00 using the concentration-mass 
relation from lMaccio et al.l |200^ for a WMAP5 cosmology. 



'Kauffmann et al."'2004; Tanaka et al. 2004:'Weinman n et al 
J006 ; Blanton c t al. 2005: Hansen et al. 2009; S kibba et al 
I2OO9I ; iBamford et aklboot) . First, we confirm that the mor- 
phological classifications described in ^2.1.11 follow the ex- 
pected morphology-density relation. Fig. [3] shows the distri- 
bution of the five morphological types as a function of E5. 
The trends with density are in agreement with the mor- 
phological trend s from other studies (e.g. iDressleJ I198GI : 
iGoto et al.l |2003| ; iBamford et al.l 120091'). Our r esults are in 
close agreement with those from Goto et aP (2003) who 
find that early and intermediate disc galaxies (Sa and SO) 
dominate in almost all environments, while late-type discs 
drop off quickly at high densities. In our sample, classical 
bulge galaxies dominate at all but the lowest densities. The 
fraction of classical bulge host galaxies rises quickly above 
E5 ~ 1 Mpc~^. Figure [3] is missing roughly 1300 galaxies 
which have axis ratios below 0.2. These galaxies are ex- 
cluded from the bulge-fdisc catalogue in L12 because our 
exponential disc model is unsuitable for edge-on discs. In- 
cluding these galaxies slightly decreases the fraction of el- 
liptical galaxies, but the changes are small and the overall 

trends remain the same. 

Note that although both iGoto et all (l2003l l and this 
work use E5 to measure local environment, the numerical 
values cannot be directly compared since E5 depends on the 
lower flux limit of the galax y catalogue used a nd the width 
of the redshift interval used. iGoto et al.l (|2003l l measure the 
distance to the fifth nearest neighbour within a redshift slice 
of ±1000 km s~^, while the redshift slice we use depends on 
the galaxy's host group velocity dispersion and is typically 
smaller. 
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Figure 4. The distribution of morphological types as a function of local density for three different stellar mass bins. The different types 
are as in Fig.|3l Although the trends with morphology at constant mass are the same as those in Fig.|3l galaxy mass is strongly-correlated 
with morphology at all local densities. 



Fig. ID shows the morphology trends for galaxies di- 
vided into three stellar mass bins. The stellar mass is com- 
puted using the most likely model for each galaxy (al- 
though the difference in total stellar mass between two rea- 
sonable model fits is always small). The trends with en- 
vironment are the same for all three mass bins, but the 
fraction of early/late type galajcies is a strong function 
of mass in all environments, in ag r eement with previous 
resul t s (e.g. [K auffma nn et al.l l2004l : iBlanton fc Moustakad 
I2OO9I : iBamford et al, 20091 ). Classical bulges dominate over 
all other types of galaxies in the two highest mass bins. 

Since colour and morphology are generally correlated, 
the morphology-density relation implies a relation between 
local density and galaxy colour such that red galaxies 
dominate in high densit y regions wh ile blue galaxies are 
dominant in the field j Gomez et al.l [2 003: Hansen c t al 
20091 : iBalogh et all l2004l : iBaldrv et al.l I2OO6I : ISkibba et al 
2OO9I ). There is evidence that the correlation between 
colour and density is more fundamenta l than the corre- 
lation between mo rphology and density ('Kauffma nn et al] 
I2OO4I : IChristlein fc Zabludoff 2005; Skibba ct al. 20Qi)- We 
show the fraction of red {u — r > 2.22) galaxies as a 
function of environment and mass in Fig. [5l In agree- 
ment with previous work, we find the fraction of red galax- 
ies is an increasing function of galaxy mass and environ- 
ment. For all masses, the fraction of red galaxies increases 
sharply near S 5 ~ 1 Mpc~^, in agreement with results from 
iGomez et all (|2003l '). who find a break in the star formation 
rate-density relation at approximately the same local den- 
sity. Figs. [3] and [5] show that both the colour-density and 
morphology- density relations hold for galaxies of con stant 
stellar mass (jBamford et al.ll2009l : ISkibba et al.ll2009l ). Ad- 
ditionally, comparing the two figures demonstrates that the 
red fraction of galaxies increases sharply at the same local 
density as the classical bulge host and elliptical galaxy frac- 
tions. We note ellipticals are not the majority of red galaxies, 
but only contribute 25 per cent by number; the remainder 
of red galaxies are classical bulge hosts, with a small contri- 
bution (5 per cent) from unclassifiable galaxies. 

In the following section, we explore the colour-density 
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Figure 5. The fraction of red galaxies (u — r > 2.22) as a function 
of S5 for galaxies in three mass bins. The fraction of red galaxies is 
highly dependent on galaxy mass, but turns up sharply at S5 ^ 1 
Mpc~2. 



and morphology-density relations for classical bulge host 
galaxies. The colours of these galaxies follow the same trends 
as the colours of the galaxy population as a whole. Us- 
ing bulge-|-disc decompositions, we can study the colour- 
density relations for bulges and discs separately, and deter- 
mine whether changes in galaxy morphology (e.g. B/T ra- 
tio) or changes in the component colours drive the increase 
in the fraction of galaxies with red integrated colour as a 
function of E5. 



4 BULGE AND DISC PROPERTIES WITH E5 

Many studies show that a galaxy's p roperties are largely 
determined by a its stellar m ass (e.g. iBrinchmann fc Ellid 
'2000!; 'Kauffmann et al' '2004") an d dark m atter host halo 
mass (e.g. Blanton et al . 2006.: ,Weinmann et al,.2006 ). Since 
the most massive haloes are strongly clustered, galaxy mass 
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and environment are strongly correlated. However, there 
are correlations with environment at fixe d stellar mafis (e.g . 
Balogh et al] 12004 iBamford et al] l2009l: [Peiig et all I2OIII: 



Griitzbauch et aimoill : ICooper et all I2OI0I : Icibinel et all 
2012I ). and these are the trends we explore. Below, we show 
trends with environment at fixed bulge mass; however, our 
results are essentially unchanged if we choose to hold total 
galaxy stellar mass fixed. In this section, we focus on our 
sample of 12500 galaxies which have both a classical bulge 
and a disc. When calculating medians and Spearman rank 
correlation coefficients, we weight each galaxy by its proba- 
bility of having a classical bulge and a disc, as explained in 

The left panel in Fig. |6] shows the colours of galaxies 
with classical bulges and discs as a function of E5. We di- 
vide the galaxies into quartiles based on bulge mass. We re- 
move galaxies with the most and least massive 0.5 per cent 
of bulges in order to eliminate galaxies with large errors in 
modelled colours and masses. The mass bin divisions are 
at logMbuige/M© = [9.3,10.1,10.3,10.6,11.2]. These mass 
bins are then divided into six bins of equal galaxy num- 
ber at different E5 . Thus, each point in Fig. [6] represents 
1/24 of the sample (~ 520 galaxies). These points are the 
weighted median colour in each bin. We define the weighted 
median as the value for which the sum of weights for val- 
ues smaller than the median is equal to the sum of weights 
for values larger than the median. The thin lines in Fig. [S] 
denote the similarly-defined weighted inter-quartile ranges. 
Although there is significant scatter in the colour, the trends 
with environment are statistically significant The trend in 
total galaxy colour is simply the colour-density relation for 
classical bulge-|-disc galaxies. The weighted Spearman rank 
coefficients are given in Table [21 We also calculate the best- 
fitting linear slopes of the relation (see Table [S]). In order 
to reduce the errors on the slope, we remove points offset 
from the fitted relation by more than 3 standard deviations 
and then re-fit the relation using the slightly smaller sample. 
Despite being statistically significant, the change in median 
integrated colour is small, ~ 0.01 — 0.03 per dex in E5. The 
change is largest for galaxies which host lo w mass bulges. 
This i s in agreement with the conclusions of lBamford et all 
(I2OO9'), who show that the colour-density relation is most 
significant for lower mass galaxies. 

The middle and right panels of Fig. [6] show the colour 
density relation for the bulge and disc components of 
these galaxies separately. It is immediately evident that the 
change in disc colour with E5 must contribute significantly 
to the change in total galaxy colour. In fact, we demonstrate 
below that the change in disc colour is the only contribution 
to the change in total colour. In the following subsections, 
we will discuss the changes in bulge properties ( H4.ip . and 
the changes in disc properties f ^4.2p as functions of E5. 



4.1 Bulges 

The classical bulges in Fig. [6] were selected to have 
large 4000A breaks, similar to ellipticals. Elliptical galaxy 



^ By statistically significant, we mean that the probability that 
the abscissa and ordinate are uncorrelatcd is less than 0.3 per 
cent, corresponding to a 3cr result for Gaussian distributions. 



Table 2. Spearman rank correlation coefficients (ps) and linear 
slopes for galaxies shown in Fig. |6] The p-value for ps is the 
probability that the colours and S5 are uncorrelated. We only 
report values larger than 10~^. The errors in the slope are the la 
errors. 



log Mh,,,]„^/Mp, 


PS 


p-value, PS 


slope 




Total g 


— r vs. log S5 




9.3-10.1 


0.17 


0.0 


0.024 ± 0.001 


10.1-10.3 


0.14 


0.0 


0.013 ± 0.001 


10.3-10.6 


0.13 


0.0 


0.012 ±0.001 


10.6-11.2 


0.11 


0.0 


0.010 ±0.001 




Bulge g 


— r vs. log S5 




9.3-10.1 


-0.04 


3.7 X 10-2 


-0.006 ± 0.002 


10.1-10.3 


-0.10 


1.1 X 10"^ 


-0.012 ±0.002 


10.3-10.6 


-0.07 


2.6 X 10-* 


-0.006 ± 0.002 


10.6-11.2 


-0.07 


2.4 X 10-* 


-0.007 ±0.002 




Disc g - 


- r vs. log E5 




9.3-10.1 


0.18 


0.0 


0.022 ± 0.002 


10.1-10.3 


0.13 


0.0 


0.017 ±0.002 


10.3-10.6 


0.12 


0.0 


0.019 ±0.002 


10.6-11.2 


0.09 


9.1 X 10-^ 


0.019 ±0.002 


Bulgc-less 


0.11 


0.0 


0.017 ±0.001 



colours depend on galaxy mass (e.g. Ide Vaucouleurd [l96ll : 
iBrinchmann &c Ellij 2000l). but only weakly, if at all, on en- 



vironment (e.g. [Dressier et al.l 19871; Bernardi et al.l l2003bl : 
iBalogh et al]|2004l : iTrager et af]|2008l : [Hansen et a"l]|2009l ). 

The central panel in Fig. [6] shows that the median bulge 
colours lie on the red sequence (g — r ~ 0.8), and that there 
is a small increase (A(gr — r) ~ 0.03) in median colour from 
the lowest mass bulges to the highest mass bulges. The weak- 
ness of this trend may be due in part to large scatter in the 
model bulge colours compared to the scatter in model ellip- 
tical colours (see fig. 32 in L12). 

In addition to trends in colour with bulge mass, there 
is a weak but statistically significant anti-correlation be- 
tween environment and bulge colour; bulges in low den- 
sity regions are redder than those in higher density regions. 
The Spearman rank correlation coefficients for this trend 
are list ed in Table [H This result is in contrast to the re- 
sults of [Hudson eTail (|2010l ) who find no variation in me- 
dian bulge colour with increasing local density. This trend 
is strongest for the second lowest bulge mass bin (red cir- 
cles). Since pseudo-bulges are typically blue and are not as 
strongly clustered as classical bulge galaxies (see Figs. [31 and 
[4} , contamination from pseudo-bulges will tend to fiatten the 
trend in bulge colour for the lowest mass bulges, possibly ex- 
plaining the weaker trend for the lowest mass bulges. 

Although the bulge colours decrease as a function of 
environment, there is no significant correlation between E5 
and SDSS fibre colours, and the correlation between E5 and 
Dn(4000) is in the opposite sense (see Fig. [SJ. However, if 
we r eplace the L12 bulge +disc decompositions with those 
from [Simard et al.[ ([201 1[ ) , the same trend in bulge colour 
with density is recovered. 

We propose that the decrease in bulge colour with envi- 
ronment is an artefact of 2-dimensional bulge+disc decom- 
positions. Specifically, the colours of small bulges in blue 
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Figure 6. Colours of bulge+disc galaxies as a function of environment. The galaxies are divided into four bins of equal number based on 
bulge stellar mass. The trends in total galaxy g — r, bulge g — r, and disc g — r are shown in the three panels. The points and thick lines 
indicate the median colours in each bulge mass bin, while the same colour thin lines indicate the inter-quartile ranges. In the disc (right) 
panel, the grey dashed lines and crosses show the median and inter-quartile range in disc colour for bulge-less (Mbuigc ~ 0) galaxies, 
including pseudo-bulge galaxies and unclassified galaxies. 
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Figure 7. Colours of classical bulges as a function of S5 using 
a model which suppresses the disc flux in the central regions. 
These models essentially eliminate the anti-correlation between 
bulge colour and E5 seen in the middle panel of Fig.|6]The points 
and thick lines indicate the median colour in each bin, while the 
same colour thin lines indicate the inter-quartile ranges. 



discs are biased redwards by the bulge-|-disc decomposition. 
The disc model used in L12 (and Smiard_ et aL (2011)) is 
an exponential profile, which is centrally peaked. Therefore, 
the disc can contribute significantly to the central flux of 
the galaxy. Moreover, since the disc is usually bluer than 
the bulge (especially in low density regions), the bulge-fdisc 
decomposition will attribute blue light from the central re- 
gion to the disc, thus making the bulge appear redder. This 
effect will be largest for the smallest bulges in the bluest 
discs. Since the bulge-to-total flux ratio and disc colour are 
both increasing functions of environment, the bias in bulge 
colour will be largest for small bulges in low density envi- 
ronments. 



We test this explanation by refitting classical bulge 
galaxies with models that suppress the disc flux in the cen- 
tral region. In the new models, the disc flux goes rapidly to 
zero within one bulge scale radius, /disc oc (r/r-buigc)/[l + 
('■/'"buige)^]"'^^- Keeping all the other model parameters 
fixed, we linearly scale the bulge and disc flux of this new 
model to fit each galaxy. By design, the new fits have larger 
bulge-to-total flux ratios, and the differences are largest for 
galaxies with small bulges. The resulting trends in bulge 
colour with environment for these disc-suppressed models 
are shown in Fig. [7] There is only a statistically significant 
trend with mass for galaxies with 10.2 < log M, ^ 10.5 
(note that because the bulge and disc masses have been re- 
computed, the mass quartiles in Fig. [7] are different from 
those in Fig. |6|. The choice of parametrisation for the cen- 
tral region of the disc plays an important role in the bulge 
colours. 



In general, it is unknown whether the stellar disc con- 
tinues unchanged through the bulge or if the disc only exists 
outside the central bulge. In the Sombrero galaxy, photom- 
etry and spectroscopy show an inner cutoff for the stellar 
disc, suggesting the disc-suppressed models may be a bette r 
choice for early-type disc galaxies (|Emsellem et al.l 1 19961 ) . 
In this work, we will continue to use the centrally-peaked 
model for the disc, knowing that the colours of small bulges 
will be biased redwards. 



This bias will also occur in the opposite case, when a 
small disc surrounds a large, red bulge. However, we do not 
expect the bias in disc colours of galaxies with large bulges 
to be as severe. Since the majority of the bulge flux comes 
from regions above and below the disc, there is no physical 
basis for the bias in disc colour, as in the case for small 
bulges. 
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Figure 8. Top: The 4000A break measured by SDSS as a function 
of S5. The sample includes both classical bulges and ellipticals. 
The different colours/symbols denote different bulge masses, and 
are the same as in Fig. [B] The points denote medians and the 
thin lines represent the inter-quartile ranges. All trends are sta- 
tistically significant. Note that the different colours/symbols do 
not all represent the same number of galaxies (as in Fig. (gJ since 
we include ellipticals in these plots. Bottom: Same as abo ve for 
galaxy stellar metallicity computed bv lGallazzi et aP l|2005l V This 
sample is based on SDSS DR4 and, therefore, only contains half 
the galaxies plotted in the top panel. 



4.1.1 Spectroscopic Properties 

Recent work has shown statistically significant positive 
correlations between stellar age and loc a l density (e.g. 
Thomas et al." 2005; Bcrnardi ct al. 2001 [ci emens et al.l 
2006. : .Smith et al. 2006; Copper ct al. 2_010|) and between 



metallicity and local density (jCooper et al. 201G| ) for early 
type, or red sequence, galaxies. Our sample of classical 
bulge+disc galaxies contains both early type, passively 
evolving, galaxies and galaxies with ongoing star formation 
in their discs. However, we expect all classical bulges and 
ellipticals to follow the same relations with stellar mass and 
environment. 

In order to examine the stellar population properties of 
classical bulges more closely, we use the SDSS fibre spec- 
tra and the line indices and metallicities reported in the 
MPA/JHU SDSS spectroscopic catalogue. At low redshift, 
the fibre spectra are dominated by bulge stellar light, but 
the mixture of bulge and disc light in the fibre is a func- 
tion of redshift. These aperture effects may influence corre- 
lations between measured line indices and E5. For classical 
bulge+disc galaxies in our sample, the median {B/T)r in 
the inner 3 arcsec is 0.74, and eighty per cent of the galaxies 
have {B/T)r within 3 arcsec larger than 0.5. For galaxies 
at 2: > 0.04, the median (3 arcsec) {B/T)r is only slightly 
smaller at 0.72. Therefore, the 3 arcsec SDSS spectroscopic 
fibres are dominated by bulge stellar light, even at the high- 
est redshift {z — 0.05) in our sample. 

Since our sample is in the regime where angular size 
scales linearly with distance, the spectroscopic measure- 



ments will include more disc flux at higher redshift. This 
effect is partially countered by a slight bias toward larger 
bulges at higher redshifts; in this sample, the average phys- 
ical bulge size increases by 8 per cent from z = 0.02 to 
z — 0.05, due to the fact that small bulges at high red- 
shift are more difficult to accurately fit. However, there is no 
statistically significant trend in E5 with redshift. Therefore, 
even though the the spectroscopic properties do change with 
redshift, the correlations between spectroscopic properties 
and E5 will be unaffected, since S5 and z are uncorrelated. 

The to p panel of Fig. [8] shows the 4000A break 
(D„(4000)) jsalogh et al.|[T999 ') as a function of Eg for four 
different bulge masses. In this figure, we include both clas- 
sical bulges and ellipticals, since we expect the stellar pop- 
ulations of cl assical bulges and ellipt icals to be the same at 
a given mass iMac Arthur et al.ll201(]| '). Indeed, excluding el- 
lipticals from Fig. [8] does not noticeably alter the results. 
The four mass bins used are the same as in Fig. [S] The in- 
crease in D„(4000) as a function of density is statistically 
significant for all bulge masses, but the trend is strongest 
for the lowest mass bulges. These results imply that bulges 
in high density environments have had less recent star for- 
mation than those in lower density environments. The same 
results are obtained if we plot the line index H5a as a proxy 
for stellar age (the equivalent width of H(5 is anti-correlated 
with environment) (|Kauffmann et al]|2003l ). The increase in 
stellar age as a function of local density at fixed bulge mass 
is in agreement with previous results (iTrager et al.l [2OO0I: 



Kauffman n et al 
2006,; Bernardi et al 



2004; Thomas ct al. 2005; Cl emens et al 
2006 ; Smith et al, 2006 ; Cooper et al 



20ld . but see lThomas et al.ll2007l '). ICooper et all (|20ld ) note 
that the trend in age of red sequence galaxies with densit y 
is evidence for galaxy assembly bias (jCroton et al.l |2007| ) ; 
namely, older galaxies are more strongly clustered than 
younger ones at fixed stellar mass. 

The lower panel of Fig. [8] shows stellar metallicity as 
a function of S5. The values for \ogZ/ZQ are taken from 
[Callazzi ct al. (2005, 20060 The metallicities are only com- 
puted for galaxies from SDSS data release 4, which in- 
cludes approximatel y half of the sample used in this paper. 
iGallazzi et ahl (|2005l ) compute th e stellar ages and metallici - 
ties by fitting model spectra from'Bruz ual fc CharlotI (|2003t ) 
to a combination of iron and magnesium line indices, the 
4OOOA break, and three Balmer lines. The trend in stellar 
metallicity at fixed bulge mass is statistically significant for 
all but the second lowest mass bin (red circles). The trends 
are strongest for the highest mass bin, where the metallicity 
changes by 0.01 dex per decade in E5. Despite being statis- 
tically significant, the change in metallicity is quite small. 

The increasing trend in stellar ages and metallicities is 
orthogonal to the ag e- metallicity degeneracy (|Worthevll994l : 
IGallazzi et al.l l2005V which makes the results more robust. 
The small increase in metalli city is in agreement with re- 
suhs fromlCooper et all l|2010l ') (but see lThomas et al.ll2005l : 
ISmith et eiL 20061 ). as well as studies of the gas-phase metal- 
licity which show an increase in metal licity in star-formin g 
galaxies as a function of local density jCooper et al.l [20081 ). 
These results demonstrate that bulges in high density re- 



* http:/ /www. mpa-garching.mpg.de/SDSS/DR4/ 
Data/ stellarmet . html 
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gion s formed earlier an d with higher star formation rates 
(see ICooper et all I2OI0I ) . The increases in stellar age and 
metallicity with increasing density suggest that classical 
bulges should be redder in higher density regions. This dis- 
agrees with the results from where we find that bulge 
colours may even be slightly bluer in high density regions. 
This disagreement is likely due to biases in our model bulge 
colours. 



4.2 Discs 

In the previous section, we show that the change in bulge 
colours with local density is small, and any statistically sig- 
nificant change is probably due to our choice of disc profile 
for galaxies with small bulges. On the other hand, the trends 
in disc colour with density are statistically significant. These 
are shown in the right panel of Fig. [6] The Spearman rank 
coefficients and slopes of the linear relation between log E5 
and {g — r)diac are given in Table [21 Although the correla- 
tions are statistically significant, the changes in disc colour 
are small; the colour only increases by ^ 0.015 mag per dex 
in E5. However, this change is 2 — 3 times larger than the 
change in colour measured for bulges. For all four bulge mass 
bins, the slope of the relation between log E5 and disc colour 
is comparable or larger than the slope in total galaxy colour, 
which suggests the change in integrated galaxy colour can be 
fully explained by the change in disc colour. We will return 
to this conclusion in H4.2.1I 

As with the bulge colours, there is significant scatter in 
disc colour at fixed bulge mass; the inter-quartile range of 
{g — r)disc colours is typically 0.1 — 0.2 mag. This scatter is 
due to the scatter in disc mass at fixed bulge mass and fixed 
S5. Disc mas s and disc colour are strongly correlated (e.g. 
Ide Jong|[l993 ). so a range of disc masses will yield a range of 
disc colours. We show below that the correlation with disc 
mass and environment does not contribute to the trend in 
disc colour with environment. In Fig. [6l we plot the rela- 
tion between E5 and disc colour for bulge-less discs, includ- 
ing pseudo-bulges and unclassifiable galaxies (grey crosses 
and dashed lines). This relation is offset to lower values of 
E5 since bulge-less discs are typically found in lower den- 
sity environments. None the less, discs in galaxies without 
prominent bulges follow the same colour-density relation as 
discs with large classical bulges. 

Like the change in total galaxy colour, the change in 
disc colour is largest for discs around the smallest bulges. 
This does not extend to bulge-less disc galaxies, for which 
the relation between log E5 and {g — r)disc is not as steep 
as the relation for galaxies with low, but non-zero, mass 
bulges (see Fig. [6] and Table [2|. However, since each bulge 
mass bin includes a large range of disc masses, this does not 
contradict the observation that the colour- density relation 
is strongest for low- (t otal) mass galaxies (jBamford et al.l 
l2009l : lTasca et al.ll2009l ). 

The trend of increasing disc colour with environment 
is a signature of star formation being halted in denser en- 
vironments. Our results are in agreement with earlier stud- 
ies of discs becoming redder in cluster environments (e.g. 
Hashimoto et al. 1 1 19981 : iMcIntosh etHI 12004 [Hudson et al.l 



els for removing gas from galaxies in dense environments 
include ram-pressure stripping by the intra-group medium 
(|Gunn &: Gott|[l972l ). and the removal of the hot halo gas 
supply around disc galaxies (strangulation) ([Larson et al.l 
IT98O). These two mechanisms have different timescales for 
shutting off star formation; ram-pressure stripping almost 
immediately ends star formation, while galaxies losing their 
halo gas rese rvoirs undergo an exponential decay in star for - 
mation rate ([Balogh et al.[[2000l :l van den Bosch et al.|[2008l ). 
From the disc fading shown here, these two processes are in- 
distinguishable. We return to the differences between ram- 
pressure stripping and strangulation in S[S|when discussing 
the correlation between disc colour and E5 for rich and poor 
groups, separately. 

Fig. [5| shows that the most massive bulges in our sam- 
ple also have the bluest median disc colours (black crosses 
in Fig. [6|. This seemingly contradicts observations that find 
redder , and presumably more massive, bulges have redder 
discs (de Jong 'l994'; 'Peletie r fc Balcellj [l996[ : [Wvse et al.1 
1997; Cameron ct al. 2009). However, Fig. [6l and [7] show 
little change in bulge colour as a function of bulge mass. 
Furthermore, while we separate the sample based on bulge 
mass, the disc mass does not monotonically increase with 
bulge mass. The trend of decreasing disc colour with in- 
creasing bulge mass exists even if the inclination correction 
is removed (see Fig. I12[) . This trend may be partially ex- 
plained by the same modelling artefact which affects bulge 
colours (34T). In this case, the subtraction of a large, red 
bulge from an image will leave behind an abnormally blue 
disc. This modelling effect will be strongest for gala^des with 
large red bulges, and relatively small discs. 



4.2.1 Disc mass and size 



2010) and the rise of red disc galaxies (anae mic spirals) 
with i ncreasing local den s ity (e.g. [i5resslerT98(]| : [Goto et al.[ 
[2OO3I : [Gomez et all [20031 : [Bamford et aiT[2009l ). Two mod- 



In ^4.21 we show that disc colour and environment are sig- 
nificantly correlated, thus disc fading is a significant contri- 
bution to the colour-density relation. However, morpholog- 
ical transformation could also play a role. If discs are being 
stripped at the same time they are fading, we expect that the 
mass contribution from the disc to be a decreasing function 
of E5. Since bulges are typically redder than discs, decreases 
in disc mass will lead to increases in total galaxy colour. In 
Fig. [9l we show the disc-to-total stellar mass ratio {D/T) 
for galaxies in four bulge mass bins. The mass bins are the 
same as in Fig. [Bl We find no statistically significant trend 
in D/T with E5. For the highest bulge mass bin (black xs), 
the median D/T is above 10 per cent, the cutoff D/T below 
which we classify galaxies as ellipticals. Therefore, it is un- 
likely that the constant D/T at high bulge mass is due to 
the minimum detectable D/T. Although we have separated 
the galaxies by bulge mass, changes in median bulge mass 
within each bin could affect the relation between D/T and 
E5. We confirm there is no statistically significant change in 
median bulge mass in each of the four mass bins as a func- 
tion of E5. Therefore, at fixed bulge mass, the disc mass is 
independent of local density, and the changes in total galaxy 
colour as a function of E5 are due solely to the changes in 
disc colour. 

Fig.dSlstrengthens the argument against morphological 
transformations by demonstrating there are no statistically 
significant correlations between disc half-light radius (mea- 
sured in the r band) and E5 at fixed bulge mass. Neither the 
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Figure 9. The stellar mass disc-to-total ratio for classical 
bulge-|-disc galaxies, at fixed bulge mass, as a function of S5. 
The bulge mass bins (indicated by colour and symbol) are the 
same as in Fig. [6] The thick lines and point denote the medians, 
and the thin lines denote the inter-quartile ranges for each mass 
bin. There are no statistically significant trends in disc mass at 
fixed bulge mass as a function of local density. 
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Figure 10. Same as Fig.|9] but for the disc half-light radius, -Reft . 
There are no statistically significant trends in disc i?(,ff s-t fixed 
bulge mass as a function of local density. 

disc size nor mass changes significantly with increased den- 
sity, as would be expected if galaxy harassment or mergers 
play a major role in galaxy evolution in high density envi- 
ronments. Thes e resul ts are in agreement with results from 
iMcIntosh et al.l (|2004l ). They find large differences in star 
formation rates between cluster galaxies and field galaxies, 
as evidenced by differences in disc colours and disc struc- 
tures (e.g. spiral arms) between cluster and field galaxies, 
but they do not find changes in the bulge-total-ratio distri- 
butions for field and cluster galaxies at fixed galaxy colour. 

It is important to keep in mind that Figs. l^land llUl are 
based entirely on galajcies with classical bulges and discs. If 
we include ellipticals {D/T = 0) in the sample, the median 
D/T is a very weakly decreasing function of S5, but only 
for the highest mass bulges. This is expected since ellipticals 
are more strongly clustered than typical bulge-|-disc galax- 
ies (see Fig. 2)). Similarly, if we include bulge- less galaxies in 
Fig. [9] (and separate galaxies into bins of constant total stel- 
lar mass instead of constant bulge mass), the median D/T 
is a decreasing function of E5 . By focusing on galaxies with 
a bulge and a disc, we are investigating whether there are 
morphological transformations within this population that 



help explain the transition from disc-dominated galaxies in 
the field to bulge-dominated galaxies in large groups and 
clusters. The lack of evolution in D/T with increasing S5 in- 
dicates that either higher density environments do not lead 
to disc destruction or the transition from a disc-dominated 
galaxy to a bulge-dominated galaxy occurs on very short 
timescales, and we do not observe any galaxies in this tran- 
sitional phase. 

4-2.2 Inclination Angle 

The previous section demonstrates that the changes in total 
galaxy colour for bulge-|-disc galaxies are due to changes in 
disc colour. However, these results could be biased by the in- 
clusion of ellipticals in high density regions. In this section, 
we address this concern by dividing the sample into three 
bins of disc inclination angle. If we mistakenly model the 
elliptical outskirts and halo as a disc, we would expect no 
colour change with environment for face-on discs, where con- 
tamination from ellipticals is most significant. Conversely, if 
contamination from ellipticals is zero and our sample does 
not suffer from any biases based on disc inclination, we ex- 
pect the same trends in bulge and disc properties with en- 
vironment, irrespective of disc inclination. 

Figures ITT] and [T^ show the disc colour and disc-to-total 
stellar mass ratio as a function of S5 for three bins in disc 
inclination angle. The galaxies plotted are the same as in 
previous figures (classical bulge-fdisc galaxies). The bins in 
inclination are chosen such that each panel shows the same 
number of galaxies. Assuming a flat disc, the inclination 
angle limits for each panel are (76°, 63°, 46°, 20°), where 
0° is face-on. From Fig. IA2I it is clear that galaxies in the 
leftmost panels of Figs. [TT] and 1121 with qa < 0.46, have a 
very low probability of being ellipticals, while the galaxies in 
the rightmost panel are the most likely to include ellipticals. 
The galaxies are divided into four bulge mass bins, using the 
same divisions as in Fig. [6] Note that the bulge mass bins do 
not all have the same number of galaxies, as they did above; 
at low qd (highly inclined galaxies), there are more massive 
bulges, while at high qd (face-on galaxies), there are more 
low mass bulges. However, the number of galaxies in each 
bin differs by at most ~ 200, and each bin has 800 galaxies 
on average. Therefore, any bias introduced by changes in 
average bulge mass with disc inclination will be small. 

Figure fTTl shows no statistically significant trend in D/T 
as a function of density. Furthermore, D/T is independent 
of disc inclination angle. This supports our previous conclu- 
sions. Namely, at constant bulge mass, disc mass is not a 
function of environment. 

The dependence of disc colour on environment and disc 
inclination is less straightforward. Although the rank cor- 
relation coefficients are always positive, the trends are only 
statistically significant for the low mass bulges or highly in- 
clined discs (see Table [31). This hints at contamination from 
elliptical galaxies in the highest qd and highest bulge mass 
bins. For the most inclined galaxies, where contamination 
from ellipticals is minimal, the slopes in disc colour reported 
in Table[3]agree with those reported for disc colours in Table 
[2I Additionally, for the lowest mass bulges (blue squares), 
the slope in {g — r)disc with E5 is essentially independent of 
qd- Although the trends for the higher mass bulges are not 
always statistically significant, the slopes measured are also 
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Figure 11. Disc-to-total mass ratio as a function of S5 for classical bulge-|-disc galaxies for three bins in disc inclination angle (reported 
as disc axis ratio, q^). The rightmost panel contains face-on galaxies. The different colours and symbols represent the same bulge mass 
bins as in Fig. |6] Each plot shows the same number of galaxies, 3240. The symbols denote the weighted medians, while the thin lines 
represent the inter-quartile ranges. 
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Figure 12. Same as Fig. 1111 but for the disc colour. The disc colours are not corrected for inclination. The correlation coefficients and 
fitted slopes are reported in Table |3] 



in reasonable agreement with those for the whole sample. 
Together with the lack of correlation between D/T and E5 
at fixed qd, the trends in disc colour with E5 at fixed qd rein- 
force our conclusion that the changes in bulge-|-disc galaxy 
colour, while small, are entirely due to changes in disc colour 
with E5. 

The colours plotted in Fig. [12] are not corrected for disc 
inclination, which accounts for the differences in the me- 
dian disc colour across the three panels. Without the in- 
trinsic inclination correction, the redder colours for inclined 
discs are expected. However, using the uncorrected colours 
will allow us to verify the disc inclination correction used 
in the previous sections. Fig. [T^] shows that the change in 
median colour as a function of qd is largest for the the two 
middle mass bins (red circles and green triangles). In this 



mass range, the median disc colour decreases by ~ 0.04 
mag from the highest inclination to the lowest inclination 
galaxies. For the high and low mass bins, the median disc 
colours only decrease 0.02 mag. This demonstrates that the 
intrinsic inclination correction should be a function of bulge 
mass; galaxies with very low mass bulges and very high mass 
bulges seem to suffer le ss intrinsic extinction in their discs. 
Following Mai ler et al. 

I iioog), our extinction correction ad- 
dresses the low bulge mass effect. The correction includes a 
term that depends on isT— band magnitude such that higher 
mass galaxies have a larger extinction correction. Thus, our 
extinction correction is typically too large for galaxies with 
massive bulges. Finally, the lack of extinction due to "discs" 
around high mass bulges may be due to contamination from 
ellipticals, where we expect little extinction. 
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Table 3. Spearman rank correlation coefficients (ps) and linear 



slopes for galaxies 


shown 


in Fig. \T2\ 




log iWbulge/™© 


PS 


p-value, PS 


slope 


Disc g — 


r vs. lo 


5 Eg, 0.25 ^ < 0.46 


9.3-10.1 


0.18 


0.0 


0.021 ±0.003 


10.1-10.3 


0.12 


1.4 X lO-"' 


0.014 ±0.004 


10.3-10.6 


0.14 


1.4 X 10"'^ 


0.025 ± 0.005 


1 O 119 


19 
U. IZ 


1.1 X lO-"' 


0.022 ±0.005 


Disc g — 


r vs. lo 


5 Eg, 0.46 ^ (?d < 0.69 


9.3-10.1 


0.23 


0.0 


0.028 ± 0.003 


10.1-10.3 


0.20 


0.0 


0.027 ±0.004 


10.3-10.6 


0.09 


5.0 X 10^3 


0.015 ± 0.005 


10.6-11.2 


0.07 


5.2 X 10^2 


0.012 ±0.005 


Disc g — 


r vs. lo 


5 Eg, Om^Qd 


< 0.94 


9.3-10.1 


0.16 


0.0 


0.020 ± 0.003 


10.1-10.3 


0.06 


0.10 


0.010 ± 0.004 


10.3-10.6 


0.06 


0.10 


0.009 ± 0.004 


10.6-11.2 


0.09 


5.6 X 10-3 


0.019 ± 0.004 



Reddening of discs is both a function of disc inclination 
and a function of the amount of dust in the disc. Above, we 
argue that discs in galaxies with different bulge masses suf- 
fer from different amounts of extinction. We can also use the 
differences in disc colour across the three panels in Fig. [12] 
to explore changes in the amount of extinction as a function 
of local density. In principle, the change in the correlation 
between the uncorrected disc colour and Eg as a function 
of qd could be used to measure the change in extinction 
as a function of local density. However, this requires a de- 
tailed understanding of how extinction depends on both the 
amount of dust and the disc inclination. Reddening is a com- 
plic ated, possibly non -monotonic, function of these variables 
fsee iTuff s ct al. 2004'. for one parametrisation) . Therefore, 
we limit ourselves to a simple test; if the extinction in high 
density environments is negligible, then the disc colour in the 
highest density environments should be independent of disc 
inclination. Fig. [T^ shows this is not true, with the possible 
exception of the highest bulge mass galaxies (black crosses) . 
As explained above, this mass bin is likely to contain mis- 
represented ellipticals, which would counter the change in 
disc colour with disc inclination. Therefore, any contamina- 
tion from ellipticals strengthens the argument that discs in 
high density regions are not free from extinction. 



5 OTHER ENVIRONMENT MEASURES 



As discussed in i|2.3l we compute environment metrics in 
addition to the local density measure. Eg. Below, we use 
the group richness, A'gai, as our environment measure. We 
define richness as the number of galaxies in a group above 
the group catalogue magnitude limit, A/o.i^ — —19.77. We 
do not correct A'^gai for contamination, nor do we attempt to 
include lower luminosity galaxies in A'gai. For rich groups, 
TVgai is closely related to a galaxy's host dark matter halo 
mass, and is notably different from Eg, which measures the 
local (intra-group) density around a galaxy. Since many of 
the effects of environment depend, at least indirectly, on the 
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Figure 14. Disc colours as a function of Eg. The sample is di- 
vided by host group richness. Isolated galaxies and galaxies in 
poor groups (A'gai ^ 20) are shown on the left, while galaxies in 
rich groups (A'gai > 20) are shown on the right. The bulge mass 
bins (indicated by colour and symbol) are the same as in Fig. 1131 
Each point in the left (right) plot represents ~ 730 (150) galaxies, 
There is no statistically significant trend in disc colour at fixed 
bulge mass for the larger groups and clusters (right panel). 



group potential, correlations between bulge and disc prop- 
erties and Agai may help determine which environmental ef- 
fects are most relevant. In Section [5. 2 1 we use crossing time, 
Across as a proxy for environment within relaxed, relatively 
rich groups. Unlike A'gai, tcross is different for each galajcy in 
a group. It measures how long a galaxy has been affected by 
its host group. 

In Fig. [13] we plot the trends in total, bulge, and disc 
colour as a function of group richness, Agai, instead of Eg. 
Unsurprisingly, the trends in galaxy colours with host group 
richness are similar to the trends with Eg ; total galaxy colour 
and disc colour are increasing functions of richness, while 
the bulge colour decreases slightly with increasing Agai. The 
latter effect is explained by the modelling bias we discuss in 
Additionally, we find no significant correlation between 
D/T and Agai, in agreement with the lack of correlation in 
Fig. [5] However, the plots in Fig.[T3]do indicate that galaxy 
colours (and disc colours) do not redden above Agai ~ 20. 
This is in agreement with results from previo u s studies (e.g. 



iBaloeh et al]l2004l : Ivan den Bosch et al.ll2008l ) . iBalogh et all 
( 200J) show that the fraction of red galaxies is independent 



of cluster velocity dispersion for a > 250 km , which cor- 
responds to Agai « 15 — 25. Thus, the majority of the colour 
evolution of bulge+disc galaxies occurs in smaller groups. 

We can test this hypothesis by examining the colour- 
density relation in relatively rich groups and poor groups 
separately. This test is shown in Fig. 1141 the left panel shows 
the relation between disc colour and Eg for groups with at 
most 20 members, while the right panel shows the same for 
larger groups. There is no statistically significant correla- 
tion between disc colour and Eg in the large groups, but the 
trends are statistically significant for galaxies in the field 
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Figure 13. Colours of classical bulge+disc galaxies as functions of group richness (Ng^i is the number of galaxies with Mo.i^ < —19.77 
in a group). Galaxies are binned by bulge mass, as in Fig. |6] The trends in total galaxy g — r, bulge g — r, and disc g — r are shown in the 
three panels. The points and thick lines indicate the median colours, while the same colour thin lines indicate the inter-quartile ranges. 
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Figure 15. The fractions of the various morphological types as 
functions of group richness. The decrease in elliptical galaxies in 
the highest bin may be due to the inclusion unvirialised (and 
possibly unbound) systems in the FoF group catalogue. None the 
less, there is little change in the morphological fractions for groups 
with more than ~ 20 members. 



and in smaller groups. Therefore, the transformation of disc 
colour does not require rich clusters, and star formation in 
discs is effectively halted by group density environments. 
This is in agree ment with other studies of colour and en- 
vironment (e.g. Zabludoff fc MulchaevI 19981: Gomez et al.l 



I2OO3I : ICooper et al.l I2OO6I : iBlanton fc Beriiiidl 120071 ). How 
ever, Fig. [14] does show a significant offset (> 0.01) in 
disc colour between galaxies in poor groups just below 
E5 « 1 Mpc~^ and galaxies in rich groups at E5 > 1 Mpc~^. 
This density threshold is close to the density of the inflection 
point in the red fraction seen in Fig. [5] 



Figure 16. The fractions of the various morphological types as 
functions of group dark matter halo mass, M200. The constant 
elliptical galaxy fraction agrees with results from iHovle et al.l 
12012) . 



In addition to examining the changes in colour as a 
function of environment, we show the changes in morpho- 
logical types as a function of A'gai and dark matter halo 
mass in Figs. 1151 and [161 In general, the trends in morpho- 
logical type fraction are in agreement with those in Fig. [3] 
For the largest groups, there is little change in the ellip- 
tical galax y fraction a s a fu nction of halo mass, in agree- 
ment with Hovle et al.l l|2012h . who find that the fraction of 
early type galaxies is « 0.2 a nd independ e nt of halo mass 



for masses above ~ 10^^ Mq. iHovle et al.l ||2012[) use mor - 



phologies from the Galaxy Zoo project ( Lintott et al.ll201ll ). 
which are qualitative classifications, and cannot distinguish 
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between face-on SOs and elliptical galaxies. As such, their 
definition of early type galaxies only corresponds approxi- 
mately to our definition of ellipticals. There are certainly 
some early type galaxies in our classical bulge-|-disc cate- 
gory. In Fig. 1151 the downturn in the elliptical fraction in 
rich groups is partially due to the inclusion of large, but 
not necessarily bound structures, in the FoF group cata- 
logue. We work to eliminate these non-virialised groupings 
the next section. 



5.1 Round groups 

Although we show that the change in disc colour oc- 
curs in relatively small groups, the original morphology- 
density relation was constructed using cluster - centri c dis- 
tance as a proxy for environment l|Oemleij 1 19741 ). In- 
deed, many studies compare galaxies on the out skirts of 

Dressleij 



mass i ve clusters to those at the centres (e.g. 



19801: iPostman fc GelleJ Il9 84': 'Whit more fc Gilmord 



Gomez et al.l 120031: iTreu et al. .20031 : ISmith et al 



199 



200' 



Weinmann et al.l l2006l: iTrager et aLTbOOSl: [Hudson et al l 



2010D . Ahhough our sample does not extend to significant 



redshifts, there are several massive groups included in the 
group catalogue which we use to study trends in bulge 
and disc properties with local density within groups. These 
trends may help reduce the scatter in the relations with E5 
and A'gai shown in Figs. [Ul and 1131 

In order to study trends within clusters, we first limit 
our sample to groups with at least 20 members. However, 
as discussed in H2.31 the FoF grouping algorithm does not 
necessarily identify relax;ed or even bound systems; fila- 
ments and merging galaxy groups are often identified by 
the FoF algorithm as galaxy groups. Furthermore, many 
large galaxy groups do not h ave well-defined cen tres (e.g. 
IZabludoff fc MulchaevI Il998l : ISkibba et al] I2OIII ). making 
trends in galaxy properties with distance from the group 
centre almost meaningless. In order to use environment mea- 
sures such as the crossing time or distance from the group 
centre, we select a sample of groups that have a symmetric 
projected distribution of galaxies. This eliminates groups in 
the process of forming or merging, for which the distance 
from the cluster centre has little physical meaning. 

If a group is relaxed and approximately spherical, we 
expect that galaxies should have a rotationally symmetric 
distribution on the sky. We generalise the requirement of 
rotational symmetry by allowing galaxies to be evenly dis- 
tributed along elliptical contours, instead of just circular 
ones. In practise, we accomplish this by compressing the 
coordinates along the minor axis of the group, making an 
elliptical group appear spherical. Allowing groups to be el- 
liptical in projection does not accurately represent oblate, 
prolate, or triaxial groups in projection (expect from special 
viewing angles), but it does account for some of the asym- 
metries real groups. 

We test whether a group is symmetric by comparing 
the distribution of galaxies observed for a given group, to 
a distribution drawn from a rotationally symmetric radial 
profile. The fu nctional form of the radial profile, however, is 
uncertain (e.g. lAdami et al.ll200ll '). Instead of relying on an 
analytic function for the profile, we create separate compari- 
son radial profiles for each galaxy group by fitting a smooth 
curve to the binned radial distribution of galaxies in each 



group. To limit the Poisson noise, this fit is smoothed using 
a Gaussian filter with a width (cr) equal to half the root- 
mean-squared (rms) radius, i?rms, of the group. 

We then compare the actual distribution of galaxies 
to the smoothed radial distribution using a test: = 

[(iVobs., i,j - Pi,j^^gai)^/(Pi,jA^gai)], where iVobs.,i,j is the 
number of galaxies observed in a region and PijN^ai is 
the expected number of galaxies from the rotationally sym- 
metric distribution. The denominator is the Poisson uncer- 
tainty in the number of galaxies. The size of the region 
used in the sum will afi'ect the x^ value; if the region is 
too large, the two distributions will trivially agree. We use a 
square region which is O.Siirms on a side. This is large enough 
to ensure most regions in the sum have at least one galaxy, 
but small enough to distinguish between rotationally sym- 
metric and asymmetric groups. If the galaxy distribution is 
rotationally symmetric, it will be statistically indistinguish- 
able from the smoothed radial distribution. On the other 
hand, if the galaxy distribution is markedly non-symmetric 
(e.g. bimodal), the x^ value will be large. We consider a 
group to be round if the x^ value from this test has a proba- 
bility between 0.1 and 0.9. This selects 95 of the 197 groups 
with at least 20 members in the group catalogue. Below, we 
designate these groups round groups. Of these groups, 46 
are represented in the bulge-f disc catalogue; the remainder 
are at redshifts greater than 0.05. There are 720 classi- 
cal bulge-l-disc galaxies in round groups, compared to 1470 
classical bulge-|-disc galaxies in groups with A'gai > 20. 



5.2 Intra-group trends 

In order to examine the efi^ect of intra-group environment on 
bulges and discs, we define the crossing time (Across) for each 
galaxy as the distance from the galaxy to the group centre 
divided by the line-of-sight velocity dispersion of the group. 
For a single galaxy, tcross is unphysical, since the galaxy is 
typically not moving radially with the group velocity dis- 
persion. However, for a large sample of galaxies, tcross is a 
measure of how long galaxies have been affected by the group 
environment. Crossing time is anti-correlated with local den- 
sity E5, although there is significant scatter. We have not 
divided our sample into central galaxie s and satellite galax- 
ies as previous studies have done (e.g. Ivan den Bosch et al.l 
I2OO8I : ISkibba fc Shetbll2009l : ISkibball2009l : iPeng et al.ll201ll ). 
However, the number of groups is much smaller than the 
number of galaxies in this sample, so the contribution from 
central galaxies is small. 

Fig. [T7] shows the relation between galaxy colours and 
Across for classical bulge-|-disc galaxies in round groups with 
A'gai > 20. Since these constraints greatly reduce the sample 
size ('^ 720 classical bulge-|-disc galaxies), we only divide the 
galaxies into three mass bins and only show three median 
points in tcross. The largest crossing times in our sample 
are less than 40 per cent of a Hu bble time, in agreement 
with those found at higher redshift (iGriitzbauch et al.| l201ll) 
and with the timescale fo r relaxation (,Gunn fc Gottlll972l : 
iFerguson fc Sandagelll990l ) . Crossing times shorter than the 
Hubble time indicate that the galaxies are not falling into 
a group for the first time, and have probably orbited the 
group centre at least once. 

The crossing time gives an indication of how long a 
galaxy has been affected by the group environment; galax- 
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Figure 17. Total, bulge, and disc colours as a function of crossing time for classical bulge+disc galaxies in round groups with at least 20 
members (see text for explanation). Crossing time is shown units of the Hubble time. The different colours/symbols represent different 
bulge masses. The points show the medians in three bins of tcioss while the thin lines show the intcr-quartile ranges. The filled points 
are the median colours for isolated galaxies (Ng^i = 1) of the same bulge mass. 



ies with the shortest crossing times have presumably passed 
through the centre of the group most often and have ex- 
perienced the highest average local density. In addition, a 
smaU tcross implies a high mass density within the group 
(tcross oc p-^^^); galaxies in compact groups will have smaller 
crossing times. Therefore, any trends in galaxy properties 
with tcross will be sensitive to how the environmental pro- 
cesses causing the trends depend on local density and time, 
although the exact dependencies are not straightforward. In 
addition, the median distance to the group centre in our 
sample of round groups is two thirds of the group virial ra- 
dius. The galaxies in our sample are all much closer to their 
host g roup centre than 3 — 4 virial radii, the distance at 
which I Gomez et all (|2003l ) find a large change in the star 
formation rate. Samples that extend to much larger radii 
are probably needed in order to see large trends in galaxy 
colours with group-centric distance. 

The left panel in Fig. [17] shows that the total galaxy 
colour is anti-correlated with tcross; galaxies with short group 
crossing times typically have redder colours. The correlation 
is only statistically significant for the highest mass bulges. 
If we do not limit the sample to galaxies from round groups, 
the statistically significance of all the trends decreases. The 
middle panel in Fig. [17] shows no statistically significant 
trend in bulge colour with crossing time. For the smaller 
bulges, there is a weak correlation with tcross, but, as shown 
in ^4.11 these trends are likely due to uncertainties in the disc 
model. The final panel in Fig. \T7\ shows the disc colours as 
a function of tcross- In this case, none of the correlations are 
statistically significant (for the highest mass bulges, tcross 
and disc colour are anti-correlated with 2a significance). 
These weak trends in disc colour are in agreement with the 
lack of correlation between {g — r)disc and E5 for all groups 
with more than 20 members (see Fig. I14p . There are sig- 
nificant colour differences between discs in large groups and 
discs in the field. Isolated classical bulge-|-disc galaxy colours 
are denoted by solid points in Fig. 1171 and, at least for low 
mass bulges, discs in the field are significantly bluer than 
discs in groups. This agrees with our earlier conclusion that 




Figure 18. Same as Fig. 1171 but showing the disc-to-total mass 
ratio (top) and the disc scale length (bottom) as a function of 
crossing time. The top plot shows that disc mass is correlated 
with crossing time. 

most disc colour evolution occurs in smaller gr oups (see also 
Ivan den Bosch et aLllioOsI : ICooper et al.ll2006l ). 

The lack of correlation bet ween disc colour and tcross 
is in contrast with the results of [Hudson et al.l (|201Cf ). who 
find a statistically significant correlation between disc colour 
and a galaxy's distance to the group centre. This is proba bly 
due to the selection of the samples; [Hudson et all (|20ld ) do 
not divide their sample into galaxies with bulges and discs 
only. If we include bulge-less galaxies in our sample, there 
is a statistically significant correlation with disc colour and 
crossing time. 

Fig. [T5| shows the change in disc-to-total mass ratio and 
disc scale length as a function of tcross- Unlike the case of 
the relations with E5, there are statistically significant cor- 
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relations between disc mass and tcross and between disc size 
and tcross- The positive correlation between D/T and fcross 
is statistically significant for the highest and lowest mass 
bins, and significant at the 2a level for the middle mass bin. 
For the highest mass bulges, the disc-to-total mass ratio in- 
creases by 6 per cent from fcross = to tcross = 0.3/-ffo- 
There is no change in average bulge mass as a function of 
tcross, so the change in D/T is due entirely to changes in the 
disc mass. In addition, the lower panel of Fig. llBl shows a cor- 
relation between fcross and the disc Rett for the most massive 
bulges, which is significant at the 2.3cr level. Although these 
trends with taoss are weak, they are also present if we sub- 
stitute Es for tcToss (correlations become anti-correlations). 
These trends suggest that processes active in the highest 
density regions in groups either suppress disc formation or 
destroy discs around infalling galaxies. We can expand our 
sample to include ellipticals {D/T = 0), to determine if the 
trends m D/T and tcross continue until the disc is negligible. 
In this case, the correlation between D/T and tcross for the 
highest bulge mass bin becomes insignificant. Thus, at least 
for high mass bulges, the trend in D/T with fcross does not 
extend to disc-less ellipticals. 

Together, Figs. 1171 and llSl tentativelv suggest that mor- 
phological changes play some role in the colour-density rela- 
tion for galaxies in rich groups. This seems to conflict with 
the conclusions in i|4.2l which show no signiflcant trends 
in disc mass with density. Taken together, these results sug- 
gest that star formation quenching and morphological trans- 
formation are separate physical processes and that these 
transitions take place in different environments. This agrees 
with previous studies that find an increase in the number 
of passive, red discs as a funct ion of time (jPressler et al.l 
1 19971 : iMoran et all [200 7: Bun dv et al.l |2010|\ and environ- 



ment (iGoto et all 1 2003: M cintosh et al.ll2004l : ISundv et all 
I2OO6I : iBamford et al.. ,200S ). The existence of these non- 
star forming disc galaxies demonstrates that star forma- 
tion quenchi ng occurs befo re morphological transformation. 
In addition, ISkibba et al.l |2009) note that there is a weak 
correlation between red galaxy morphology and density at 
small scales, in agreement with our findings. Namely, we 
demonstrate that morphological transformation may be tak- 
ing place in rich groups and clusters, and that these trans- 
formations are not associated with star formation quenching 
since disc colours are not correlated with taoss- Galaxy ha- 
rassment, i.e. high-speed encounters between galaxies, is a 
plausible explanation for the decrease in disc mass we ob- 
serve. The number of encounters a galaxy experiences will 
increase with galaxy number density, and galaxies with short 
crossing times are exposed to the highest average densities. 
Thus galaxy harassment would yield correlations between 
tcroaa and disc mass and disc size similar to those in Fig. [TS] 



6 SUMMARY 

In this work, we examine the changes in bulge and disc prop- 
erties for a sample of 12500 galaxies as functions of local 
projected density. Since galaxy mass and environment are 
strongly correlated, we divide the sample into bins of equal 
bulge mass in order to study any residual trends in bulge 
and disc properties with local density. Using 2-dimensional 
bulge-l-disc decompositions, we are able to study the colour- 



density and morphology-density relations for bulges and 
discs, separately. Our sample consists of galaxies with a 
classical, elliptical-like bulge surrounded by a disc. Classical 
bulges are observed to have the same characteristics as ellip- 
ticals of the same mass; the only difference is the encompass- 
ing disc. By studying the properties of classical bulge hosts 
in different environments, we can deduce if local density is 
the determining factor in whether or not a classical bulge 
acquires and retains a disc. In addition, we study whether 
the population of classical bulge+disc galaxies is undergo- 
ing a transition from disc-dominated and star-forming in 
low-density regions to bulge-dominated and passive in high- 
density regions, as is expected from the morphology-density 
and colour-density relations. 

Based on our results, we can draw two conclu- 
sions about the effects of environment on discs around 
classical bulges. First, both the colour and the mass 
of these discs changes with increasing density, but 
these changes occur at different densities and, presum- 
ably, on different time-scales. This suggests that star 
formation quenching and morphological transformation 
are caused by different physical processes, in agree- 
ment with_jnajiypreyio^ environment 
(e.g. iGoto et all l2003l: iGomez et al] l2003l: [M cintosh et al] 



2004 



Kauffmann et al.'"2004': 'Christlein & Zabludofi ,2005 



[Cooper e t al. 2006; Bundv et al. 2006; van den Bos ch et al.l 
'2008-; Bamford ct al. 20091: iSkibba et al.ll2009l : iHudson et all 
i201Qy Second, although disc properties are a function of 
environment, the changes in disc colour and disc mass are 
insufficient to explain why classical bulges have discs and el- 
lipticals do not have discs. While environment clearly affects 
disc formation and evolution around classical bulges, there 
must be other parameters that determine whether or not a 
classical bulge is surrounded by a disc. 

The first piece of evidence for our conclusions comes 
from separating the colour-density relation for bulge+disc 
galaxies into relations for bulges and discs. We find that 
the correlation between total galaxy g — r and local den- 
sity is due entirely to changes in disc colour with density. 
The small change in bulge colour with increasing density 
depends on our choice of disc model, but this choice does 
not noticeably affect disc colours. At fixed bulge mass, disc 
mass and disc scale length are independent of local density. 
Thus, there is no evidence for morphological transformation 
as a function of local density, but there is evidence for star 
formation quenching in discs with increased local density. 
Morphological transformation and star formation quench- 
ing must occur separately, and any process responsible for 
star formation quenching cannot dramatically alter stellar 
discs. 

Even for low bulge masses, where the colour-density 
relation is strongest, the change in disc colour with increas- 
ing density is small; (g — r)disc increases by ~ 0.05 mag 
over two decades in projected galaxy number density. This 
colour change is not enough to explain the colour-density 
relation for the entire galaxy population, which is twice 
as steep. Therefore, the correlations between stellar mass 
and density and between stellar mass and galaxy colour 
contribute significantly to the colour-density relation for 
classical bulge+disc galaxies. We expect processes inter- 
nal to galaxies, which may depend on stellar mass, are at 
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least as important in moderating star formation as external, 
environment-related, processes. 

The colour-density relation for discs is not linear. By 
dividing the sample into groups with more or fewer than 
20 members, we show that the relation between density 
and disc colour disappears for galaxies in rich groups, re- 
gardless of bulge mass. This is in agreement with results 
that star formation is quenched well before galajcies en- 
ter massive clusters (iGomez et al.ll2003l : ICooper et al.ll20o"^ : 
iMcIntosh et al.|[2004l '). It also puts limits on the physical pro- 
cesses responsible for stopping star formation. For example, 
ram-pressure stripping requires cluster-scale velocities and 
is unlikely to account for star formation cessation in smaller 
groups. Strangulation, on the other hand, is effective at rel- 
atively low densities and is, therefore, a possible cause of 
star formation quenching in discs. Although our analysis is 
not sensitive to colour gradients in discs, changes in these 
gradients as a function of environment may put additional 
constraints on the physical processes responsible for disc fad- 
ing as a function of local density (see .Roediger ct a l. 2011). 

After determining that star formation quenching occurs 
outside rich groups, we examine trends in bulge and disc 
properties within rich groups. Instead of local projected den- 
sity, we use the group crossing time, tcross, as a measure of 
environment. Since tcross depends on distance to the group 
centre, we restrict our sample to "round" groups, which have 
a smooth angular distribution of galaxies. Even for these 
groups, the luminosity-weighted group centre may not be 
physically relevant; it is not unusual fo r the most massive 
galaxy to be offset from the group centre jSkibba et al.|[201ll : 
liCeorge et ai]l2012l '). The average crossing time for galaxies in 
this restricted sample is roughly one third of a Hubble time, 
suggesting most of these galaxies have already orbited their 
host group at least once. We find no statistically significant 
correlations with disc colour and fcross, but we do find mod- 
erately significant (2 — 3a) correlations between disc-to-total 
mass ratio (D/T) and tcross, and between disc scale length 
and f cross, at fixed bulge mass. The trends m D/T are not 
a result of correlations between group richness and tcross, or 
between richness and D/T. In addition, the correlations are 
also present at the same significance if we replace tcross with 
local projected density. Like the changes in disc colour as 
a function of projected density, the changes in disc mass as 
a function of tcross are small; the median D/T decreases by 
less than 10 per cent from the largest to the smallest crossing 
times. None the less, we infer that bulge-|-disc galaxies do 
undergo morphological transformations in large groups, but 
they do not undergo star formation quenching at the same 
time. This is direct evidence for the separation of morpho- 
logical transformation and star formation quenching. 

The above results demonstrate that environment has 
two distinct effects on the discs surrounding classical bulges. 
First, disc star formation is truncated in (relatively) poor 
groups, leading to the colour-density relation for discs 
(Fig. [6]). Although gas-stripping requires relatively high ve- 
locities, tidal interactions and heating by the intra-group 
median may remove a disc's outer halo gas supply (strangu- 
lation) , thus quench ing star formation over several gigayears 
jLarson et al.lll980l ). This preprocessing of galaxies in small 
grou ps has been su ggested before as the origin of SO galaxies 
(e.g. iDressler et al.i,199 7V The second effect is a morpholog- 
ical transformation. As these quenched galaxies enter higher 



density environments in larger groups, the stellar disc is dis- 
rupted over several orbits, leading to the observed D/T- 
icross correla tion. This disrupti on may be caused by galaxy 
harassment (|Moore et al.|[l99^ '). which is most effective for 
galaxies in high density regions with short crossing times. 

Although environment does affect discs around classi- 
cal bulges, these effects are small. Therefore, while statisti- 
cally significant, the changes in disc colour with density and 
disc mass with crossing time are insufficient to explain the 
full range of the colour-density relation and morphology- 
density relation. Furthermore, since the changes are small, 
present-day local density cannot be the determining factor 
in whether or not a classical bulge has a disc. There must 
be other processes, unrelated to present-day environment, 
which regulate disc formation around classical bulges. None 
the less, because classical bulges and ellipticals of the same 
mass seem to have the same formation history, the pro- 
cesses that regulate disc formation are probably external to 
the galaxy. By only examining trends with density at fixed 
bulge stellar mass, we cannot explore trends in bulge mass 
with density. Environmental processes, s uch as ram pressure 
stripp ing, may drive bulge growth (e.g. iTonnesen fc BrvanI 
I2OO9I I. We plan to explore the effects of environment on 
bulges at fixed total stellar mass in a later paper. 

One way to further explore disc formation around clas- 
sical bulges would be to look for evolution in discs around 
classical bulges as a function of redshift. In this paper, we 
present correlations between disc properties with environ- 
ment in the local universe. In drawing conclusions, we have 
assumed these correlations are signatures of the evolution of 
bulges and discs as galaxies move from low density environ- 
ments to high density environments. However, these conclu- 
sions need to be supplemented with observations of bulges 
and discs at higher redshifts. The comparison data needed 
for such studies is easily available from space-based optical 
and near-infrared s urveys. For e x ample , using data from the 
CANDELS survey, iBruce et akl (|2012l ) presents bulge-f-disc 
decompositions for massive galaxies beyond z = 1. Careful 
comparisons of bulges and discs at different redshifts and 
in different environments will better constrain what effect 
environment has on the evolution of galaxies and their com- 
ponent bulges and discs. 
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APPENDIX A: CLASSIFICATION OF 
GALAXIES 

Below, we describe the galaxy classification scheme briefly 
outlined in tj2.1.1l We explain the classification of each 
type of galaxy, paying special attention to the probabilistic 
separation of ellipticals, classical bulge+disc galaxies, and 
pseudo-bulge+disc galaxies. 

(i) Bulge-less galaxies: First, we identify galaxies which 
are best fit by a disc-only model. We assign them a prob- 
ability of being fit by a single exponential equal to one. 
These bulge-less galaxies are selected because they have an 
de Vaucouleurs bulge-to-total luminosity ratio in the r band 
{{B /T)r) smaller than 10 per cent, or because there is no 
statistically significant improvement in the fit by adding a 
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Figure Al. The absolute (left) and apparent (right) magnitude 
distributions for galaxies from the L12 sample classified using the 
classification scheme described in ^2.1. H and appendix|X] Galaxies 
fit with a Sersic profile are unclassifiable. They make up 17 per 
cent of the sample from L12, and less than 10 per cent of the 
bright {Mr < —19.77) subsample used in this work. 

bulge component. These galaxies account for 37 per cent 
of the L12 sample (26, 882 galaxies). These galaxies are not 
strictly bulge- less, but include galaxies with bulges too small 
to detect using bulge-fdisc decompositions. Fig. lAll shows 
that these galaxies are predominately intrinsically (and ap- 
parently) faint, as expected for disc galaxies with, at most, 
very small bulges. Because we use a bright (Afr < —19.77) 
subset of galaxies from L12, bulge- less galaxies are a smaller 
fraction (28 per cent) of our sample than of the whole L12 
sample. 

(ii) Ellipticals and red disc galaxies: Next, we select 
a sample of elliptical galaxies. As above, galaxies with 
{B/T)r > 0.9 are assumed to be di sc-less ellipt i cal g alax- 
ies. As discussed in L12 (see also lAllen et all [2006 ). 2- 
dimensional bulge-|-disc models of elliptical galaxies usually 
consist of a de Vaucouleurs profile and a low surface bright- 
ness "disc" . This model "disc" component has several possi- 
ble origins: the outer halo of ellipticals, a Sersic index larger 
than 4, and/or inadequate sky subtraction around bright 
galaxies in SDSS. This model "disc" component makes el- 
liptical galaxies indistinguishable from dustless, face-on, red 
classical bulge-|-disc galaxies based on photometry alone. 
Therefore, we must use other information to separate red 
ellipticals from face-on red classical bulge-l-disc galaxies. Be- 
low, we describe a method for separating face-on bulge-|-disc 
galaxies from ellipticals. We only apply this separation to 
galaxies with red colours for both the model de Vaucouleur s 
bulge and model disc: u — r > 2.22 (jStrateva et al.l [20011 ). 
This colour cut ensures that the outer exponential halo 
around elliptical galaxies will not be star-forming. However, 
we cannot identify blue ellipticals in our sample using this 
colour cut. 

Despite their similarities, we can separate red bulge-|-disc 
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Figure A2. Top row: The disc axis ratio distribution for red bulge+red disc galaxies and ellipticals divided into three stellar mass 
bins. These galaxies are selected to have red (u — r > 2.22) nj, = 4 bulges and discs. The dotted curve is the best-fitting disc axis ratio 
distribution for red B+D galaxies (flat above > 0.1) and the dashed curve is a Gaussian fit to the distribution of disc axis ratios from 
bulge+disc models fit to ellipticals. The sum of these curves is shown by the solid line, which is fit to the histogram. Bottom row: The 
probability that a red galaxy with a given axis ratio and mass is a red classical bulge+red disc galaxy as opposed to an elliptical. 



galaxies and ellipticals based on the statistics of the mod- 
elled disc axis ratios, qa- If the discs are randomly ori- 
ented, then qd should be uniformly distributed for val- 
ues larger than the disc scale height to scale length ratio, 
qz- The probability density function is given by f{qd) = 
Id/ \J [id ~ ~ We use q^ = 0. 1, which is smalle r 

than the average value qz = 0.14±0.04 ferege l et al]l2002[ ). 
However, any qz less than 0.2 has a negligible effect on the 
results. If the bulge-|-disc models are applied to ellipticals, 
then the distribution of qd will be skewed to higher values, 
mimicking face-on discs. Although it is impossible to distin- 
guish a dustless, face-on classical bulge-|-disc galaxy from an 
elliptical, we can assign each galaxy a probability of being 
an elliptical or bulge-l-disc galaxy based on its model qd- 

Fig. I A2I shows the statistical separation of red bulge+disc 
and ellipticals for galaxies divided into three stellar mass 
bins. The upper panels show the distribution of qd for galax- 
ies with red (u — r > 2.22) model de Vaucouleurs bulge and 
exponential disc colours, and 0.1 < {B/T)r < 0.9. Each 
distribution of axis ratios has two contributions, one from 
disc galaxies, which follows the distribution for randomly 
oriented discs, and one from ellipticals, which we model as 
a Gaussian distribution in qd with an arbitrary centroid and 
width. We separately fit each distribution of qd in Fig. IA2I 
with a linear combination of these two probability density 
functions. Since small measured axis ratios are often due to 
poor fits, we only use galaxies with qd > 0.25 for the fit- 
ting. The fractional contribution from each function gives 
the probability that a galaxy with a given qd is an ellipti- 
cal or a red classical bulge+red disc galaxy (lower panels 
Fig. IA2|) . For low mass galaxies, the distribution of axis ra- 
tios does not show any contribution from ellipticals. At any 
stellar mass, a small value for qd implies a galaxy has a real 
disc. For 0.6 < gd < 0.8, the probability that a galaxy is an 
elliptical rises to ~ 50 per cent for the highest mass galaxies. 



This method of distinguishing between bulge+disc galax- 
ies and ellipticals does not allow us to definitively assign a 
galaxy to either category, but we can still examine the prop- 
erties of the whole sample, simply by weighting each galaxy 
by its likelihood of being an red bulge+red disc galaxy or 
an elliptical. Counting all galaxies with B /T > 0.9 as el- 
lipticals, we find that 9.4 per cent (6800 galaxies) of the 
L12 sample are red bulge+red disc galaxies and 7.3 per cent 
(5300 galaxies) are ellipticals. 

When selecting galax;ies with red bulge and disc compo- 
nents for Fig. IA21 we first correct the bulge and disc colours 
for inclination in order to eliminate edge-on, dusty galaxies. 
However, if the inclination corrections are too large, there 
will be fewer galax;ies with small qd, and the fraction of ellip- 
ticals will be enhanced. Therefore, we iteratively adjust the 
inclination corrections in u and r such that, after identifying 
elliptical galaxies, there is no residual trend in corrected disc 
colour with disc inclination for red bulge+red disc galaxies. 
The final adjustments make the inclination corrections in u 
and r smaller than those in L12, but the adjustments are 
minimal. They change the number of ellipticals in sample 
by less than 1 per cent of the total sample (500 galaxies). If 
we used no inclination correction, the number of ellipticals 
in the sample would decrease by 1000 galaxies. 

(iii) Classical and pseudo-bulges: After identifying ellip- 
ticals and red classical bulge+disc galaxies, we classify the 
remaining galaxies as either classical bulge hosts or pseudo- 
bulge hosts. We separate classical bulges and pseudo-bulges 
based on the age of the central stellar populations. Fis herl 
(|2006l ) show that pseudo-bulges selected by morphology have 
higher central star formation rates than classical bulges, 
in agreement with the notion that pseudo-bulges form via 
secular processes in discs which drive gas inw ards and en- 
hanc e the central star formation (|Kormendv fc KennicuttI 
l2004h . Fig. [Ml shows the distribution of the 4000A break 
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Figure A3. Left: The distribution of the 4000A break 
(Dn(4000)) for classical (magenta) and pseudo- (cyan) bulges. 
The total distribution (black line) is fit with two Gaussians in 
order to separate the the two types of bulges. Right: The distri- 
bution of the disc axis ratios for classical and pseudo-bulges. The 
excess of edge-on discs around classical bulges is due to the in- 
clusion of edge-on red B-l-D galaxies. The lack of edge on discs 
around pseudo-bulges is due to the difficulty of detecting a fiat 
pseudo-bulge in an edge-on disc. 



strength for the classical and pseudo-bulge galaxies. These 
values are taken from the MPA/JHU catalogue of measured 
spectroscopic quantities from SDSS (|Tremonti et al.ll2004l : 
lAihara et al. 2011.fl. The value plotte d is the narrow def i- 
nition of the 4000A break, D„(4000) (|Balogh et al.lll999l ). 
The SDSS spectrograph uses 3 arcsec fibres; two-thirds of 
the galaxies shown have a bulge-to-total fiux ratio within 
3 arcsec larger than 0.5. Therefore, the fibre spectroscopic 
quantities are typically dominated by the stellar fight from 
the bulge. 

Fig. IA3I shows there are two populations of bulges, those 
with recent star formation (within ~ IGyr) and small 
D„ (4000) and those without recent star formation and large 
D„ (4000) . This suggests that D„ (4000) can be used to sep- 
arate quiescent classical bulges from star-forming pseudo- 
bulges. We fit the distribution of D„ (4000) with two Gaus- 
sians and assign each galaxy a probability of having a classi- 
cal bulge or pseudo-bulge based on the ratio of the Gaussians 
at a given Z)„(4000). Once again, an individual galaxy's type 
is unknown, but the statistical properties of the whole sam- 
ple can be examined. Sixty per cent of the bulges in Fig. lA3l 
are classical bulges. Classical bulge galaxies make up 18 per 
cent of the total L12 sample, while pseudo-bulges represent 
12 percent. 

The right panel of Fig. IA3I shows the distribution of disc 
axis ratios for classical and pseudo-bulge host discs. This 
is one test of our separation of pseudo-bulges and classi- 
cal bulges, as we expect the distributions to match the fiat 
distribution for randomly oriented discs. Although the dis- 
tribution of Qd for classical bulge host discs together with 



pseudo-bulge host discs is flat, the separate distributions for 
classical and pseudo-bulge hosts are not. The pseudo-bulge 
host discs have axis ratios which are too large. This is ex- 
pected, since pseudo-bulges are typically flattened, making 
them difficult to detect in edge-on discs. Edge-on pseudo 
bulge hosts are more likely to be considered bulge-less galax- 
ies. 

The slight excess of low axis ratio discs around classical 
bulges is due to the inclusion of inclined, but dust-poor, 
red bulge-|-red disc galaxies. When selecting red bulge+disc 
galax;ies and ellipticals above, we included an inclination cor- 
rection, although we expect that some bulge+disc (and most 
ellipticals) will be essentially dust-free. Correcting these 
galaxies for inclination shifts their colours bluewards, re- 
moving them from the red sample. This is especially true 
for highly- inclined (low qd) galaxies, where the colour cor- 
rections are largest. Therefore, the excess of small qd classi- 
cal bulges galaxies consists mainly of red classical bulge-|-red 
disc galaxies. Since all these galaxies host classical bulges, 
it is not important to accurately separate them based on 
colour. 

The physical origins of classical and pseudo-bulges are 
very different. Classical bulges are thought to have formed 
by the same mechanisms as ellipticals, while pseudo-bulges 
arise fro m secular processes in discs (dis c instabilities, bars, 
etc) (e.g. iKormendv fc Kennicuttl |2004| ) . Furthermore, the 
absolute magnitude distribution of pseudo-bulge galaxies in 
Fig. lAll is more similar to that of bulge-less and unclassifi- 
able galaxies than it is to classical bulge hosts. In this work, 
we will consider galaxies with pseudo-bulges as a subset of 
bulge-less disc galaxies. Because pseudo-bulges galaxies are 
relatively faint, they are a small fraction (8 per cent) of 
the bright {M,. < —19.77) subsample from L12 used be- 
low. Therefore, the exclusion of pseudo-bulge hosts from 
bulge-|-disc galaxies does significantly affect our results. 

(iv) Unclassifiable galaxies: The unclassifiable galaxies in 
our sample are modelled by a single Sersic profile. These 
are galaxies which are not well fit by any of the other mod- 
els, i.e. they have model bulges which are larger than their 
discs, or they have fluxes in g, r, ot i consistent with zero 
for the bulge or disc (or both). They are given a probability 
of being flt by a Sersic profile equal to unity. Of the 72658 
galaxies in L12, 12459 (17 per cent) are deemed unclassi- 
fiable. The average unclassifiable galaxy is 0.4 magnitudes 
fainter than the average galaxy in the sample (see Fig. IA1|) . 
Seventy-five percent of unclassifiable galaxies have a Sersic 
index less than 2.3 and the same fraction lie in the blue cloud 
{u — r < 2.22). These galaxies are probably disc-like irreg- 
ulars, which are unlikely to have a well-defined bulge and 
disc. The remaining 25 per cent of unclassifiable galaxies 
are mostly merger remnants, starbursts, and other compli- 
cated morphologies. None the less, because the majority of 
unclassifiable galaxies exhibit disc-like properties, we group 
them with other bulge-less galaxies described above. Since 
the sample used here is a bright subsample of the L12 sam- 
ple, the unclassifiable galaxies make up a smaller fraction 
(less than 10 per cent) of the bright subsample than of the 
whole L12 sample. 
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